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1.  INTRODUCTION  AND  CHAPTER  OVERVIEW 

This  chapter  considers  basic  research  related  to  the  extreme  environment 
of  an  aircraft  engine  and  the  use  of  Thermal  Bander  Coatings  (TBC’s) 
to  ameliorate  the  effects  of  extreme  temperature  cycling  on  metal  engine 
components.  The  failure  of  these  TBC’s  is  a  serious  technological  problem; 
one  that,  if  solved,  should  greatly  increase  the  fuel  efficiency  and  operating 
lifetimes  of  airplane  engines.  These  TBC’s  arc  comprised  of  ceramics, 
with  favorably  low  thermal  conductivity,  deposited  on  the  engine  metals. 
Accordingly,  we  arc  concerned  with  the  characterization  of  Metal-Ceramic 
(M/C)  interfaces  at  a  fundamental  level.  In  this  chapter,  we  attempt  to  provide 
an  overview  of  experimental  techniques  for  characterizing  M/C  interfaces. 
However,  since  we  arc  theorists,  much  of  the  review  is  focused  on  providing 
a  detailed,  critical  analysis  of  theoretical  methods  in  use  today  to  study  such 
systems.  We  also  give  examples  from  our  own  modeling  at  the  atomic  level 
that  has  yielded  some  insights  into  the  interfacial  behavior  of  TBC’s. 

Let  us  outline  the  problem  caused  by  the  extreme  combustion  environ¬ 
ment.  The  ideal  engine  would  operate  at  very  high  temperatures  without 
failure  in  order  to  have  the  highest  fuel  efficiency.  Typically,  the  combustion 
gas  is  held  at  temperatures  above  1370  °C,  while  the  metal  superalloys  that 
constitute  the  engine  components  have  melting  points  ranging  from  1230- 
1315  °C!  This  makes  it  imperative  to  either  cool  the  metal  components  (e.g., 
by  drilling  holes  and  flowing  cool  air)  or  to  provide  thermal  protection  from 
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the  combustion  gas.1  Ideas  for  optimizing  the  cooling  techniques  and  engine 
metal  alloy  compositions  reached  a  point  of  diminishing  returns  at  least  10 
years  ago.2  Hence,  engineers  looked  to  ceramic  materials  as  a  means  of 
providing  a  thermal  barrier  coating  that  will:  (i)  extend  the  life  of  gas  tur¬ 
bine  components,  (ii)  reduce  cooling  requirements  (thereby  decreasing  fuel 
consumption),  or  (iii)  allow  for  an  increase  in  gas  inlet  tempratures  (thereby 
increasing  thrust).3 

Ceramics  are  thermal  insulators  that  can  provide  the  thermal  barrier  de¬ 
sired.  The  challenge  of  working  with  such  ceramics  is  that  typically  they 
have  completely  different  thermochemical  properties  from  the  material  to 
which  they  are  expected  to  adhere,  namely  a  metallic  alloy.  As  a  result, 
the  usual  thermal  cycling  that  such  M-C  interfaces  undergo  tends  to  stress 
these  interfaces  to  the  point  of  fracture  and  spallation  (chipping-oif).1  It  has 
been  suggested  that  research  is  needed  to  try  to  connect  these  macroscopic 
phenomena  to  microscopic  properties,  in  order  to  make  progress  in  under¬ 
standing  how  to  design  the  best  M/C  junction.4-5  In  particular,  the  exact 
nature  of  the  interface  at  the  atomic  level  is  poorly  characterized;  the  ex¬ 
act  mechanisms  have  yet  to  be  identified  by  which  the  interface  is  formed, 
stressed,  fractured,  and  spalled.  Furthermore,  it  remains  unclear  what  roles 
oxidation  and  temperature  play  in  stabilizing  or  destabilizing  the  interface. 
Gaining  an  atomic-level  understanding  of  these  mechanisms,  should  help 
elucidate  ways  to  optimize  the  M/C  couple:  simultaneously  maximizing 
thermal  insulation  and  minimizing  spallation. 

This  chapter  initially  explores  the  TBC  from  an  experimental  point  of 
view.  In  Section  2  of  this  chapter,  we  present  a  detailed  description  of  the 
TBC  structure  and  composition.  We  explore  the  typical  chemical  make-up 
of  a  TBC,  the  formation  of  the  thermally  grown  oxide,  and  the  physical 
methods  used  in  TBC  fabrication.  Section  3  lists  a  number  of  the  experi¬ 
mental  techniques  used  to  characterize  functional  TBC’s  and  explains  some 
experimental  characterization  of  ideal  interfaces.  Section  4  investigates  the 
problem  of  TBC  spallation.  This  includes  identifying  likely  perpetrators 
of  spallation,  describing  proposed  spallation  mechanisms,  looking  at  some 
features  that  complicate  spallation  studies,  and  expounding  on  the  issue  of 
adhesion.  We  also  mention  some  of  the  techniques  used  for  measuring 
stress,  fracture,  and  spallation  at  the  end  of  this  section.  After  looking  at 
TBC’s  from  an  experimental  standpoint,  we  turn  our  discussion  to  present 
applicable  theory.  We  describe  atomistic  modeling  approaches  in  Section  6 
and  how  those  can  be  applied  to  M/C  interfaces.  This  includes  exploring  the 
use  of  cluster  and  slab  models,  the  idea  of  interface  stoichiometry,  and  the 
problem  of  lattice  misfit.  We  then  consider  the  available  ab  initio  techniques 
for  modeling  M/C  interfaces.  This  is  broken  down  to  more  specific  meth- 
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ods  including  quantum  chemical  approaches,  density  functional  theory,  the 
approximation  of  the  Harris  functional,  and  tight  binding  schemes.  We  also 
present  results  obtained  from  a  number  of  density  functional  theory  studies, 
including  our  own,  as  well  as  tight  binding  predictions.  Section  7  and  8 
review  theory  designed  to  handle  larger  systems  than  arc  computationally 
feasible  with  the  ab  initio  techniques.  Finally,  we  offer  brief  conclusions 
that  attempt  to  look  to  the  future  for  ways  of  enhancing  understanding  of 
and  optimization  of  thermal  barrier  coatings. 

2.  CREATING  A  TBC  SYSTEM 
2.1.  Commonly  Used  Materials 

When  creating  a  TBC,  generally  a  top  coat  and  a  bond  coat  layer  must  be 
deposited.  The  top  coat  serves  as  the  insulator  and  the  bond  coat  mediates 
contact  between  the  top  coat  and  metal  alloy  substrate.  The  nickel-based 
“superalloy”  substrate  is  a  real  pot  pourri  of  elements,  consisting  of  Ni,  Co, 
Cr,  Mo,  Al,  Ta,  Ti,  C,  Zr,  and  B  (where  the  non-Ni  elements  are  present  at 
the  few  to  hundredths  of  a  weight  percent  level.)6  Yttria-Stabilized  Zirco- 
nia  (YSZ)  is  a  favored  top  coat  material.  Pure  zirconia  has  relatively  high 
strength,  wear  resistance,  and  fracture  toughness.  Likewise,  it  exhibits  an 
extremely  low  thermal  conductivity.  In  fact,  excluding  Pyrex  glass,  the  ther¬ 
mal  conductivity  of  zirconia  is  lower  than  any  other  engineering  ceramic  by 
over  an  order  of  magnitude.6  In  practical  terms,  this  property  allows  even  a 
thin  zirconia  film  (less  than  one  millimeter  thickness)  to  potentially  reduce 
the  temperature  of  the  underlying  alloy  several  hundred  degrees  Celsius.7 
Furthermore,  the  linear  thermal  expansion  coefficient  and  elastic  modulus  of 
zirconia  (especially  the  tetragonal  phase)  are  well-matched  to  several  popular 
nickel-based  superalloys,  compared  to  possible  alternative  ceramics.  These 
properties  arc  essential  to  coating  survival  during  thermal  cycling.7-9  Unfor¬ 
tunately,  zirconia  also  exhibits  polymorphism  as  a  function  of  temperature, 
primarily  between  the  monoclinic,  tetragonal,  and  cubic  phases.  Although 
the  polymorphism  can  be  exploited  to  inhibit  crack  propagation  via  volume 
changes  that  occur  upon  transformation,71011  unregulated  polymorphism 
can  be  detrimental  to  TBC  functionality.  As  a  result,  other  cubic  oxides 
arc  added  to  control  or  eliminate  polymorphism.  Partial  stabilization  of  the 
tetragonal  phase  across  the  relevant  temperature  range  affords  some  phase 
control  yet  preserves  the  valuable  inhibition  of  stress-induced  micro-crack 
propagation.  Greater  than  8.5  mole  percent  Y2O3  dopant  produces  fully- 
stabilized  cubic  zirconia  while  a  2-8.5  percent  Y2O3  concentration  creates 
a  partially  stabilized  (tetragonal)  form;12  similarly,  adding  Ce02,  CaO,  or 
MgO  can  generate  stabilized  zirconia.13-16  Conversely,  using  UO2  as  the 
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stabilizing  oxide  is  less  effective  than  the  other  oxides  studied. 1718  However, 
Ti02  added  to  Y2O3  or  CeC>2  stabilized  Z1O2  has  been  shown  to  produce  zir- 
conia  polycrystals  with  favorable  properties.19  20  Although  a  percentage  of 
the  partially  stabilized Z1O2  tends  to  remain  in  the  monoclinic  phase  through¬ 
out  thermal  cycling,  the  initially  formed  cubic  phase  converts  to  tetragonal 
after  thermal  treatment  in  air.21-22  Despite  the  favorable  properties  of  a  YSZ 
top  coat,  the  difference  in  thermomechanical  properties  between  a  YSZ  top 
coat  and  a  metal-alloy  substrate  is  enough  to  require  the  introduction  of  an 
intermediate  layer.  This  bond  coat  is  important  for  adhesion  and  grading  the 
thermal  expansion  mismatch  between  the  top  coat  and  substrate.  A  typical 
bond  coat  contains  nickel,  chromium,  aluminum,  and  yttrium,  with  nickel 
as  the  primary  element  for  nickel-alloy  substrate  applications.  With  other 
alloys,  it  is  sometimes  desirable  to  use  iron  or  cobalt  in  place  of  the  nickel 
in  the  bond  coat.21-24 

2.2.  Formation  of  the  Thermally  Grown  Oxide 

A  third  layer  present  in  a  TBC  is  the  Thermally  Grown  Oxide  (TGO). 
Bond  coat  oxidation  can  reduce  TBC  adhesion  to  the  substrate.  In  most 
cases,  the  oxidation  products  begin  to  form  even  prior  to  top  coat  deposition. 
It  is  standard  practice,  in  zirconia  film  creation  chambers,  to  supply  oxygen 
to  ensure  stoichiometry  of  the  zirconia  layer.  Moreover,  the  bond  coat 
continues  to  be  oxidized  once  the  zirconia  layer  is  in  place,  since  zirconia 
readily  conducts  oxygen  ions  at  high  temperature.25  Accordingly,  the  TGO 
layer  between  the  bond  coat  and  top  coat  thickens  with  thermal  cycling. 
Although  this  layer  is  generally  thin  compared  to  the  TBC  layer,  it  can  lead 
to  large  stresses  in  the  system  due  to  significant  thermal  expansion  mismatch 
of  the  TGO  and  bond  coat. 

The  growth  of  the  TGO  falls  into  two  primary  regimes.  Fast  initial 
growth  of  non-protective  oxides  is  followed  by  continuous  protective  scale 
growth,  largely  controlled  by  diffusion.  The  primary  protective  scale  formed 
with  bond  coat  oxidation  is  alumina,  AI2O3.  For  a  bond  coat  composed  of 
NiCrAlY,  the  oxidation  also  produces  Ni(Cr,Al)204  spinels,  Y2O3,  NiO, 
AIYO3,  and  AI5Y3O12  bands  oriented  perpendicular  to  the  bond  coat.23 
In  the  initial  oxidation  stages,  X-ray  diffraction  of  a  FeCrAlY  bond  coat 
indicate  that  both  FcCi^CH  and  C12O3  form  in  that  system.24 

Oxide  scale  studies  are  complicated  further  by  the  coexistence  of  various 
oxide  phases.  For  instance,  at  early  times,  both  0-alumina,  a  transient  phase 
formed  when  y- alumina  is  heated,  and  a-alumina  arc  present.  Amorphous 
alumina  transforms  through  a  series  of  metastable  phases  until  the  more 
stable  a-alumina  dominates  at  long  oxidation  times.26  Although  thicknesses 
of  the  layers  vary,  some  fairly  typical  values  arc  250-500  pm  for  the  top 
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coat,  100-150  urn  for  the  bond  coat  (this  varies  since  some  bond  coats  are 
designed  for  surface  roughness  and  with  several  gradation  layers),  and  about 
an  order  of  magnitude  less  than  these  for  the  thermally  grown  oxide.23  Fig.  1 
displays  a  schematic  cross-section  of  a  TBC. 


Figure  1  Cross-section  of  a  TBC  (not  pictured  to  scale) 


2.3.  Depositing  the  TBC  and  Bond  Coat 

Several  methods  effectively  apply  thin  coatings  to  substrates.  One  process 
commonly  employed  for  TBC  fabrication  is  thermal  plasma  spraying.  For 
this  technique,  the  bond  coat  is  deposited  with  low  pressure  plasma  spraying 
and  the  Zr02  top  coat  is  then  created  using  atmospheric  plasma  spraying.27 
A  potential  problem  with  this  method  is  coating  fracture  due  to  weak  in- 
terlamellar  adhesion,  resulting  from  discontinuities  within  the  solidifying 
process.28  Accordingly,  plasma  sprayed  coatings  tend  to  fail  within  the  TBC 
itself.  Alternatively,  electron-beam  physical-vapor  deposition  can  be  used 
for  TBC  fabrication.  One  drawback  of  this  method  is  possible  spallation,  via 
separation  at  the  bond  coat  interface  during  cooling.26  Nevertheless,  success¬ 
ful  coatings  of  this  type  possess  columnar  micro-structure,  supplying  some 
of  the  stress  compliance  required  by  the  TBC.  Furthermore,  this  method  is 
suited  for  creating  a  graded  AI2O3-Z1O2  layer.  It  is  hoped  that  a  gradual 
transition  from  the  bond  coat  oxide  to  the  top  coat  would  improve  the  life 
of  the  TBC;29  however,  more  work  is  needed  to  understand  the  practical 
effectiveness  of  such  a  layer.30  Sputter  deposition  is  another  useful  means 
for  producing  thin  coatings.31’32  Each  deposition  method  has  its  own  set  of 
pros  and  cons,  hence  this  variety  in  fabrication  techniques  continues. 
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3.  CHARACTERIZING  THE  TBC 

3.1.  Characterization  Techniques 

Advanced  characterization  methods  are  required  to  gain  a  microscopic- 
level  understanding  of  TBC’s.  High-Resolution  Transmission  Electron 
Microscopy  (HRTEM)  and  Scanning  Electron  Microscopy  (SEM)  are  fre¬ 
quently  used  to  explore  the  structure  of  the  coating  surface  and  the  metal- 
oxide  interface.22, 33,34  SEM  also  provides  a  qualitative  thickness  measure 
of  oxide  growth.  The  coating’s  porosity  can  be  measured  by  image  analysis 
techniques  and  mercury  porosimetry.  Use  of  the  latter  can  pose  difficulties 
since  large  pores  arc  tilled  with  mercury  prior  to  applying  pressure.  Con¬ 
versely,  image  analysis  is  ineffective  for  measuring  very  small  pores  and 
microcracks.35  X-ray  diffraction  provides  a  means  to  determine  the  oxide 
species  and  phases.22  Likewise,  field-emission  SEM  and  energy  dispersive 
spectroscopy  aid  in  characterizing  oxidation  products.23  Microindentation 
tests  determine  the  elastic  modulus  of  the  prepared  coatings;  however,  due  to 
inherent  inhomogeneity  in  the  coatings,  results  from  this  method  must  be  ap¬ 
propriately  averaged.35  Thermal  conductivity  measures  have  been  provided 
by  laser  flash  methods  that  measure  diffusivity  and  specific  heat  measure¬ 
ments  from  differential  scanning  calorimetry.36,37  Alternatively,  diffusivity 
values  can  be  obtained  with  a  multiproperty  apparatus  that  measures  tem¬ 
perature  gradients,  sample  geometry,  and  heat  flux.38  The  single-wavelength 
pyrometer  method  gathers  spectral  effective  emissivities.  Although  thermo¬ 
couples  could  also  be  used,  the  single-wavelength  pyrometer  technique  has 
the  superior  attributes  of  not  requiring  surface  contact  during  temperature 
measurements,  immunity  to  electromagnetic  interference  from  the  surround¬ 
ings,  and  no  system  perturbation  during  measurements.39  The  use  of  Auger 
Electron  Spectroscopy  (AES),  X-ray  photoelectron  spectroscopy  (XPS),  and 
Auger  parameter  (a’)  analysis  provide  means  to  explore  the  nickel-alumina 
interface  formation.  Furthermore,  these  techniques  arc  suited  for  determin¬ 
ing  ionicity  and  growth  mechanisms  at  metal-oxide  interfaces.40 

3.2.  Characterization  of  Ideal  Interfaces 

At  high  temperatures,  interdiffusion  leads  to  the  formation  of  new  phases 
at  the  M/C  interface.4  For  example,  Qin  and  Derby41  used  optical  and  elec¬ 
tron  microscopy  to  characterize  the  strength  of  Ni/Zr02  and  Ni/NiO/ZrCU 
interfaces.  They  found  that  the  strongest  interfaces  were  formed  by  anneal¬ 
ing  Ni/Zr02  in  air,  which  allowed  a  thin  layer  of  NiO  to  form  that  helped 
adhesion  to  the  ceramic.  Qin  and  Derby  also  studied  the  formation  of  an 
interface  between  Z1O2  and  a  Ni(Cr)  alloy,  where  a  mixed  oxide  with  a 
spinel  structure  was  formed  at  the  interface.42  The  reaction  at  the  interface 
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appeared  to  be  accompanied  by  local  melting  after  interdiffusion,  although 
it  is  not  clear  which  elements  were  actually  diffusing.  Wagner  et  al.43  ex¬ 
amined  charge  flow  in  a  NilZrC^IZr  cell  and  showed  that  at  1273  K  with  no 
applied  voltage,  oxygen  ions  move  from  the  Ni  to  the  Zr  electrode,  forming 
Ni  -jZr  and  N17Z12  at  the  Ni  electrode  while  reducing  Zr02  to  Zr.  At  the  Zr 
electrode,  monoclinic  Z1O2  forms.  Thus,  oxygen  ion  diffusion  appeal's  to 
be  an  important  process  by  which  adhesion  takes  place. 

As  mentioned  earlier,  AI2O3  (TGO)  also  forms  at  the  TBC-bond  coat- 
superalloy  junction.  Trumble  and  Riihle44  showed  by  Transmission  Electron 
Microscopy  (TEM)  that  at  1 390  0  C,  pure  Ni  does  not  form  a  spinel  (MAI2  O4) 
at  the  interface  between  Ni  and  AI2O3.  However,  even  0.07  percent  oxygen 
in  the  Ni  will  induce  formation  of  a  thin  spinel  layer  at  the  interface  (with 
no  NiO  intermediate  required),  where  the  kinetics  appeal'  to  be  controlled  by 
oxygen  diffusion.  Shear  strength  measurements  by  Loh  et  al.45  determined 
that,  under  conditions  where  the  spinel  NiAl204  is  formed,  fracture  occurs 
along  the  spinel-Ni  interface.  It  is  still  controversial  as  to  whether  formation 
of  the  spinel  compound  at  these  interfaces  actually  helps  or  hinders  adhesion. 
Zhong  and  Ohuchi,46  using  X-ray  photoelectron  spectroscopy,  and  later 
Brydson  et  al.,5  using  spatially  resolved  transmission  Electron  Energy  Loss 
Spectroscopy  (EELS)  and  High  Resolution  Electron  Microscopy  (HREM), 
determined  that  the  interface  between  Ni  and  AI2  O3  forms  direct  Ni- Al  bonds 
under  reducing  conditions  at  high  temperatures.  They  presented  evidence 
for  a  M3AI  phase  at  the  interface,  which  they  suggested  was  formed  by 
Al  diffusion  into  the  Ni,  and  that  this  phase  provides  a  driving  force  for 
the  formation  of  Ni-Al  bonds  rather  than  Ni-0  bonds.  They  suggested  that 
formation  of  any  spinel  phase  containing  Ni-0  bonds  needs  to  be  minimized 
because  the  spinel  is  brittle,  an  exactly  opposite  conclusion  to  that  reached 
by  Qin  and  Derby  for  the  Ni-Zr02  interface. 

In  addition  to  M/C  interactions,  characterization  of  Zr02-Al203  is  inter¬ 
esting  as  well,  because  of  the  AI2O3  that  forms  between  the  bond  coat  and 
the  TBC.  Aita  and  coworkers  have  grown  and  characterized  with  HREM 
nanolaminate  films  that  alternate  between  polycrystalline  Z1O2  and  AI2O3 
layers.  They  have  shown  that  Z1O2  grows  with  the  close-packed  tetrag- 
onal(lll)  or  monoclinic(lll)  surfaces  parallel  to  the  substrate.47  Interest¬ 
ingly,  we  recently  calculated  those  surfaces  to  be  the  most  stable  for  each  of 
these  two  phases  48  They  also  observed  that  the  amount  of  t-Zr02  increases 
with  decreasing  Z1O2  film  thickness  and  that  the  crystallites  grow  in  this 
tetragonal  phase  up  to  a  critical  thickness  of  about  6.0  nm,  at  which  point 
additional  Z1O2  converts  the  crystallite  into  the  monoclinic  phase.49  In  other 
work,  they  showed  that  the  stress-induced  transformation  of  tetragonal  to 
monoclinic  zirconia  is  limited  to  nanometer-scale  regions  in  these  nanolam- 
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inates.50  In  Section  6.2.2.2,  we  discuss  our  own  theoretical  rationalization  of 
these  observations.  The  role  of  alumina  seems  merely  to  confine  the  size  of 
Z1O2  crystallites  formed.  Our  own  calculations  on  Z1O2-AI2O3  (discussed 
in  Section  6.2.2.2)  are  consistent  with  these  observations. 

4.  TBC  FAILURE 

Differences  in  crystal  structure,  size,  and  thermal  expansion  coefficients 
between  the  top  coat,  bond  coat,  and  substrate  introduce  strain  into  the  TBC. 
The  thermal  expansion  coefficient  of  the  top  coat  is  generally  lower  than 
that  of  the  bond  coat  or  the  alloy  substrate.  As  mentioned  earlier,  a  properly 
designed  bond  coat  serves  as  a  graded  thermal  expansion  layer  to  reduce 
the  strain  caused  by  thermal  expansion  mismatch.  Unfortunately,  thermal 
cycling  in  air  enhances  the  mismatch  and  increases  strain  resulting  from 
interfacial  damage  during  oxidation.24  With  repeated  cycling,  this  strain 
contributes  to  coating  failure. 

4.1.  Likely  Culprits 

TBC  failure  involves  a  number  of  contributing  factors.  Obviously,  the 
strains  introduced  with  thermal  cycling  are  a  major  area  of  concern.  The 
buildup  of  the  TGO,  migration  and  segregation  of  the  components,  and 
phase  transitions  of  the  TGO  and  top  coat  may  each  contribute  to  the  coating 
failure.  Even  the  type  of  porosity  in  the  coating,  largely  dependent  on  the 
coating  procedure  used,  can  affect  the  coating  lifetime.  A  great  deal  of 
experimental  work  connected  with  TBC’s  exists  in  the  literature.  Ernst  pub¬ 
lished  an  excellent  review  that  covers  many  experiments  performed  before 
1995.4  Due  to  inherent  complexities  of  TBC’s,  studies  to  elucidate  the  fail¬ 
ure  mechanisms  generally  concentrate  on  only  certain  aspects  contributing 
to  failure.  The  subtle  relationships  of  all  the  factors  are  not  fully  understood; 
hence,  further  modeling  and  experiments  are  needed  to  gain  a  more  thorough 
understanding  of  present  weaknesses  and  possible  future  improvements. 

4.2.  Spallation  Mechanisms 

Spallation  is  the  process  by  which  the  TBC  peels  off  of  the  substrate;  and 
naturally,  after  the  coating  has  spalled,  the  continuous  thermal  protection 
layer  no  longer  exists.  Likewise,  the  spalled  fragments  can  block  gas  flow, 
contaminate  products,  and  permit  corrosion.  The  failure  of  the  TBC  leading 
to  spallation  appeal's  closely  related  to  the  damage  process  within  the  TGO 
layer.  Unfortunately,  due  to  the  highly  complex  nature  of  the  problem,  there 
is  limited  understanding  of  how  growth  kinetics  of  the  TGO  and  microcrack 
damage  affect  TBC  lifetime.  Three  primary  scenarios  have  been  proposed 
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for  the  spallation  mechanism.51  The  first  of  these  states  that  a  buckling 
effect  can  result  from  planar  compressive  stresses  within  the  ceramic  layer. 
It  seems  this  is  a  plausible  mechanism  provided  there  exists  an  interfacial 
delamination  crack,  formed  due  to  local  conditions  creating  both  out-of- 
plane  and  shear  tensions.  Recent  microscopy  data  indicate  spallation  can 
proceed  via  this  buckling  mechanism  when  the  delamination  crack  is  at 
least  sixteen  times  the  TBC  thickness.51'52  However,  for  thick  oxides,  the 
buckling  mode  is  not  viable.  Those  systems  may  undergo  a  surface  wedging 
effect.53  This  second  model  relies  on  the  development  of  a  through-thickness 
shear  crack  in  the  TBC  due  to  compression.54  Finally,  a  more  recently 
proposed  mechanism  depends  upon  a  void  formation  under  the  TBC  and 
subsequent  folding  effects  which  may  lead  to  cracking.55  This  effect  is 
also  known  as  “wrinkling”  or  “rumpling.”  Naturally,  in  addition  to  voids 
formed  through  thermal  cycling,  the  pore-types  initially  created  in  the  TBC 
may  influence  the  eventual  spallation  mode.  Although  each  model  appeal's 
plausible  in  its  description  of  the  final  spallation  effect,  an  understanding 
that  would  elucidate  dynamic  evolution  conditions  leading  to  spallation  is 
desired.  Fig.  2  displays  schematic  diagrams  of  the  three  spallation  models. 


Buckling 


Wrinkling 

Figure  2  Spallation  Mechanisms  of  the  TBC 

4.3.  Complicating  Features  of  Spallation 

Although  an  understanding  of  spallation  mechanisms  is  useful,  it  seems 
TBC  and  substrate  delamination  results  from  a  highly  complex  relationship. 
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which  includes  bond  coat  oxidation,  micro-crack  evolution,  and  progressive 
buckling.51  Fracture  surfaces  created  by  TBC  spallation  suggest  that  failure 
occurs  in  the  TBC-bond  coat  interfacial  region.56  A  piezospectroscopic  study 
shows  large  failure  regions  may  follow  the  oxide  grain  impressions  in  the 
bond  coat,  indicating  failure  where  the  TGO  has  grown  into  the  bond  coat 
grain  boundaries.26  Certain  morphological  instabilities  of  the  oxide -bond 
coat  interface  may  cause  the  fracture  that  joins  the  failure  regions.  These 
instabilities  could  result  in  local  normal  forces  across  the  interface  leading 
to  interfacial  fracture.  A  similar  phenomenon  occurs  with  spontaneous  spal¬ 
lation  after  cooling  a  TBC  to  room  temperature.  Consequently,  sub-critical 
crack  growth  may  be  hastened  by  environmental  factors.  Propensity  for 
moisture  to  enhance  such  crack  growth  in  metal-alumina  interfaces  further 
supports  these  ideas.26- 57- 58  Using  thermal  plasma  spraying  for  bond  coat 
deposition  may  result  in  a  layered  bond  coat  with  irregular  thickness.  Sam¬ 
ples  possessing  this  discontinuity  in  bond  coat  thickness  failed  after  thermal 
treatment  under  an  argon  atmosphere,  indicating  that  at  least  one  factor  in 
thermal  spallation  is  increased  residual  stress  on  the  YSZ  after  thermal  treat¬ 
ment.21  Thus,  although  the  buildup  of  the  TGO  may  play  a  large  role  in 
the  eventual  TBC  failure,  an  overall  spallation  mechanism  encompasses  a 
variety  of  complicating  features.51 

4.4.  Adhesion 

Composition  of  the  substrate  and  bond  coat  affects  adhesion  of  alumina 
formed  during  bond  coat  oxidation.  Of  course,  when  the  alumina  spalls, 
the  zirconia  layer  is  not  maintained;  so  the  TBC  fails.  Migration,  segrega¬ 
tion,  and  stress  generation  arc  key  areas  of  concern  in  the  thermally  grown 
oxide  region.  Migration  of  aluminum  ions  from  bond  coat  to  metal-alloy 
substrate  occurs  due  to  a  concentration  gradient.  Likewise,  migration  from 
the  substrate  and  bond  coat  into  the  thermally  grown  oxide  region  is  also 
expected.  Actually,  it  appeal's  that  limited  yttrium  ion  migration  to  the 
alumina  improves  adhesion;25  however,  contaminants  in  the  alumina  gen¬ 
erally  lead  to  additional  stresses,  which  decrease  adhesion  during  thermal 
expansion.59  The  high  temperature  diffusion  of  bulk  yttria  ions  to  the  YSZ 
surface  destabilizes  the  top  coat.60-61  This  permits  zirconia  phase  transfor¬ 
mation  to  monoclinic62  resulting  in  an  undesirable  volume  expansion  that 
may  contribute  to  spallation.  Furthermore,  an  SEM  study  shows  bond  coat 
aluminum  can  diffuse  into  the  YSZ  layer.39  Previously,  it  has  been  suggested 
that  both  neutral  and  ionic  aluminum  diffusion  into  the  YSZ  layer  induces 
spallation.63  Accordingly,  diffusion-controlled  migration  in  the  TBC  may 
promote  harmful  phase  transitions  and  enhance  thermal  stresses. 
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Although  a  concentration  gradient  allows  migration  via  diffusion,  the  oxy¬ 
gen  chemical  potential  gradient,  arising  during  oxidation,  supplies  another 
driving  force.  Studies  show  that  this  chemical  potential  gradient  penetrates 
the  top  coat,  bond  coat  and  into  the  metal  substrate.  It  affects  oxygen- 
reactive  species  in  all  these  layers,  including  those  initially  in  the  form  of 
stable  oxides,  nitrides,  carbides  and  sulfides.64-67 

Segregation  to  the  metal-TGO  interface  also  affects  oxide  adhesion.  For 
instance,  sulfur  segregates  to  this  interface,  which  can  prove  detrimental 
to  TBC  lifetime.  The  interfacial  sulfur  increases  the  thickness  of  the  oxi¬ 
dation  layer,  decreases  adhesion  of  the  oxide  layer  to  the  metal,  increases 
transformation  of  metastable  alumina  to  the  alpha  phase,  and  enhances  pore 
formation  at  the  interface  and  within  the  oxide  layer.  Increased  interfacial 
roughening  and  void  formation  results.  The  formation  of  voids  is  prob¬ 
lematic  since  voids  act  as  stress  concentration  sites  within  the  oxidation 
layer.65-68-72  De-sulfurization  of  the  sample  to  less  than  1  ppm  may  help 
prevent  spallation.73-76  However,  it  is  possible  that  the  presence  of  other 
species  (such  as  Y  and  Zr)  in  the  bond  coat  and  substrate  may  supress  the 
harmful  effects  of  the  sulfur  and  render  de-sulfurization  unnecessary. 25,65,77 
Finally,  scale  stresses  can  result  from  isothermally  generated  growth  stress 
during  the  oxide  formation  or  from  cooling  stress  due  to  thermal  expansion 
mismatch  during  thermal  cycling.78 

4.5.  Measures  of  Stress,  Fracture  and  Spallation 

In  addition  to  experimental  means  for  initially  characterizing  TBC  sys¬ 
tems,  a  number  of  techniques  arc  useful  for  investigating  stress,  fracture,  and 
spallation.  Cr3+  piezospectroscopy  is  an  optical  method,  sensitive  to  Cr3+ 
dopants  in  the  alumina,  that  permits  study  of  oxidation  stresses  through  the 
TBC.  Prior  to  this  technique,  non-destructive  study  of  the  oxide  layer  had 
proven  difficult  due  to  the  physical  location  of  the  alumina.  Fortunately, 
zirconia  is  fairly  transparent  at  the  frequencies  of  interest.  As  a  result, 
this  piezospectroscopic  technique  is  able  to  detect  stresses  and  some  phase 
differences  within  the  oxide  layer.79  Techniques  for  determining  fracture 
energy  include  double  cantilever-beam  experiments,80-81  four-point  bending 
tests,82-84  and  wedge-loaded  peel  tests.83  Furthermore,  laser  spallation  pro¬ 
vides  a  means  to  measure  tensile  strength  of  the  metal-oxide  interface.85-87 
Pulsed  lasers,  used  as  high  shock  generators,  allow  exploration  of  the  spalla¬ 
tion  process  through  simulated  high  pressure  loading.88-89  A  method  based 
upon  Thin  Layer  Activation  (TLA)  attempts  to  directly  measure  spallation. 
TLA  relies  on  the  creation  of  radionuclides  in  the  surface  layer  after  exposure 
to  a  high-energy  beam  of  charged  particles.90,91  Loss  of  activated  material 
due  to  spallation  results  in  a  decreased  y-activity  signal.  Furthermore,  it 
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is  possible  to  collect  the  spalled  material,  allowing  mass  quantification  that 
provides  an  additional  sensitive  spallation  measure.92  Recently,  effective 
measurements  of  amplitude  and  profile  of  laser  generated  stress  pulses  have 
been  made  using  a  Doppler  velocity  interferometer  system  for  any  reflector. 
An  advantage  of  time  resolved  interferometry  is  that  it  allows  the  researcher 
to  investigate  dynamic  behavior  of  the  material  by  looking  at  damage  ef¬ 
fects  on  wave  propagation.93  Dilatometry  provides  a  quantitative  means  to 
explore  interfacial  damage  due  to  oxidation,  allowing  thermal  expansion 
experiments.24  The  strain  introduced  by  thermal  expansion  mismatch  plays 
a  major  role  in  TBC  failure,  so  techniques  for  exploring  thermal  expansion 
are  very  important. 

5.  THEORETICAL  PREDICTION  OF  INTERFACE  STABIL¬ 
ITY 

Theoretical  methods  offer  the  opportunity  to  explore  structure -property 
relationships  in  ideal  metal-ceramic  interfaces.  Ultimately,  improved  un¬ 
derstanding  of  the  causal  sequence  leading  to  a  particular  interface  structure 
and  set  of  properties  would  enable  further  optimization  of  manufacturing 
parameters.  Atomistic  modeling  constitutes  the  perfect  laboratory  in  this  re¬ 
spect.  Within  the  limits  of  the  specific  approximations  used  for  interatomic 
interactions,  physical  properties  may  be  resolved  to  arbitrary  accuracy  and 
competing  effects  may  be  separated. 

Here  we  will  be  concerned  mainly  with  the  interface  structural  stability 
and  electronic  structure.  One  of  the  most  important  factors  determining 
the  interface  stability  is  the  interface  cohesive  strength.  Other  factors  may 
be  equally  important  for  stability  of  the  interface,  depending  on  the  actual 
situation,  e.g.  corrosion  resistance,  thermal  expansion  and  elastic  matching 
of  the  metal  and  ceramic,  the  flexibility  of  the  structure  to  release  stress 
buildup  during  thermal  cycling,  stability  towards  structural  degradation  by 
unwanted  interface  segregation  of  certain  elements  and  undesirable  mixed 
phase  formation  at  the  interface.  In  this  section,  we  give  an  overview  of 
strategies  by  which  theory  can  play  an  important  role  in  helping  characterize 
such  interfaces.  A  recent  and  excellent  review  on  the  theoretical  aspects 
of  the  M/C  interface  was  given  by  Finnis94  in  1996;  we  will  therefore 
concentrate  on  later  developments  in  the  understanding  of  M/C  interfaces. 

6.  ATOMISTIC  APPROACHES 

The  discipline  of  atomistic  modeling  has  proliferated  tremendously  over 
the  past  two  decades,  due  to  the  increased  capacity  of  modern  computers. 
The  basic  trade-off  in  atomistic  modeling  is  always  between  accuracy  of 
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calculated  energies  and  forces  and,  on  the  other  hand,  size  of  the  atomic 
ensemble  meant  to  model  a  macroscopic  system.  The  larger  the  ensemble 
size,  the  smaller  is  the  influence  of  miscellaneous  finite  size  effects  on 
the  physical  properties.  However,  the  larger  the  system  size,  the  more 
approximate  the  atomistic  description  necessarily  becomes. 

Which  tradeoff  to  choose  depends  on  the  situation:  some  physical  prop¬ 
erties  arc  insensitive  to  details  in  the  interatomic  potentials.  In  such  cases, 
a  model  potential  will  provide  essentially  identical  results  to  the  exact,  as 
long  as  the  model  potential  reproduces  certain  characteristic  quantities  like, 
e.g.,  elastic  constants  and  bulk  cohesive  energies. 

Physical  properties  may  be  insensitive  to  details  of  the  interatomic  poten¬ 
tial  for  a  variety  of  reasons.  For  example  at  low  temperatures,  the  atomic 
system  may  probe  a  bounded  region  of  phase  space,  most  often  the  elas¬ 
tic  regime  around  the  equilibrium  state;  and  in  this  bounded  region,  the 
interatomic  interactions  may  be  represented  well  by  a  suitable  model  po¬ 
tential.  Some  classes  of  physical  properties,  like  critical  phenomena,  are 
intrinsically  insensitive  to  many  details  in  the  interatomic  potentials,  as  they 
are  controlled  by  collective  behavior.  For  other  physical  phenomena,  like 
melting,  effects  of  fine  details  in  the  interatomic  potentials  tend  to  average 
out. 

However,  in  other  situations,  there  will  be  a  strong  sensitivity  to  details  in 
the  interatomic  potential.  This  is  the  case  when  the  physics  and  chemistry 
is  governed  by  rare  events,  like  crossing  activation  barriers  and  the  breaking 
and  forming  of  chemical  bonds.  If  the  transition  state  is  not  well  character¬ 
ized  and  included  carefully  in  the  parameterization  data  base  for  the  model 
potential,  the  activation  barrier  is  very  likely  to  be  described  erroneously 
by  the  model  potential,  leading  to  wrong  transition  rates.  This  is  where  ab 
initio  methods  arc  needed. 

In  the  sections  below,  the  particular  aspects  associated  with  atomistic 
modeling  of  M/C  interfaces  will  be  reviewed  along  with  applicable  recent 
work  done  in  this  field. 

6.1.  Structural  Models  for  M/C  Interfaces 

Atomistic  approaches  need  a  structural  model  to  represent  the  real  system. 
In  this  section,  we  will  briefly  discuss  the  most  common  choices  of  structural 
models  along  with  advantages  and  disadvantages,  in  order  to  enable  readers 
unfamiliar  with  details  to  critically  judge  results  presented  in  the  literature. 
We  will  then  discuss  issues  involved  in  choosing  the  structure  of  the  in¬ 
terface  which  involves  lattice  stoichiometry  and  misfit.  The  most  common 
structural  models  of  M/C  interfaces  fall  into  two  main  classes:  cluster  and 
slab  representations. 
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Figure  3  Schematic  picture  of  M/C  interface  structural  models. 

A  :  Cluster  model  with  vacuum  around  it. 

B  :  Dense  M/C  interface  unit  cell,  physically  corresponding  to  a  superlattice  or  sandwich 
structure. 

C  :  Slab  model  with  vacuum  between  periodic  images  perpendicular  to  the  interface.  Physi¬ 
cally,  the  structure  corresponds  to  M/C  thin  films.  The  borders  of  the  periodically  replicated 
unit  cells  in  B  and  C  are  shown  with  bold  lines. 

6.1.1.  Cluster  Models  These  models  are  usually  constructed  by  scooping 
out  a  small  representative  region  (see  Fig.3A)  containing  5-30  atoms  from 
the  real  M/C  interface.  Cluster  models  emphasize  chemistry  as  local  by 
only  considering  a  small  “active”  region.  Unfortunately,  some  early  cluster 
studies  used  unrealistically  small  clusters,  which  were  poor  representations 
of  the  real  interface  they  attempted  to  model.95  In  such  cases,  important 
long-range  effects  are  missed. 

Sometimes  the  ceramic  side  of  the  cluster  is  embedded  in  an  array  of  the 
nominal  anion  and  cation  point  charges  in  the  proper  structure,  to  emulate 
the  real  Madelung  potential  in  the  ceramic  -  this  definitely  improves  the 
realism  of  the  cluster  model.96  100  A  Green’s  function  constructed  from  the 
perfect  host  crystal  has  also  been  used  to  embed  a  ceramic  cluster.101  An 
alternative  to  anion  and  cation  point  charges  is  embedding  the  cluster  into 
an  array  of  overlapping  anion  and  cation  pseudopotentials,102' 103  thus  trying 
to  capture  some  of  the  electron-electron  interaction  with  the  surrounding 
ceramic  medium.  A  supplement  to  this  is  mechanical  embedding  of  the 
cluster  via  force  field  interactions  with  the  surrounding  substrate  lattice,104 
so-called  molecular  mechanics.  The  latter  neglects  the  perturbations  in 
interatomic  force  constants  caused  by  metallic  adsoiption  and  is  likely  to  be 
rather  inaccurate  in  the  case  of  reconstruction  or  extensive  relaxations  in  the 
interface  region. 

Cluster  models  have  other  generic  finite  size  effects.  Interfaces  are  often 
geometrically  frustrated:  for  all  metal  atoms,  there  is  a  preferred  adsoiption 
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site  on  the  ceramic  surface.  When  one  metal  atom  occupies  that  site,  the 
neighboring  metal  atom  generally  will  not  be  able  to  occupy  an  equally 
favorable  site,  due  to  the  lattice  mismatch  between  metal  and  ceramic.  This 
frustration  effect  is  not  modeled  properly,  if  the  metal  side  of  the  cluster  is 
too  small.  Another  finite  size  effect  is  that  the  electronic  density  of  states 
oscillates  significantly  with  the  cluster  size.105  This  is  critical  near  the  Fermi 
level,  because  such  oscillations  in  the  density  of  states  available  for  bonding 
to  the  ceramic  is  likely  to  produce  predictions  that  converge  slowly  with 
cluster  size.106  This  effect  may  be  countered  by  “state  preparation”105  or  by 
embedding  the  metal  cluster  in  a  periodic  metal  slab. 107  Unfortunately,  since 
many  finite  size  effects  converge  slowly  with  the  cluster  size,  unmanage¬ 
ably  large  clusters  may  be  required  if  not  embedded  properly.  Furthermore, 
charged  clusters  arc  often  considered,  due  to  the  nonstoichiometry  of  the 
ceramic  piece  of  the  cluster.  Otherwise,  the  anions  and  cations  will  not  be 
in  the  proper  charge  state,  corresponding  to  the  overall  neutral  bulk  ceramic. 
A  back-of-the-envelope  estimate  suggests  that  the  nonphysical  polarization 
induced  by  charged  fragments  may  substantially  affect  the  chemical  predic¬ 
tions.  The  one  significant  advantage  of  cluster  models  is  that  the  power  of 
ab  initio  quantum  chemical  methods  can  be  applied  and  therefore  systemat¬ 
ically  converged  results  independent  of  experiment  can  be  obtained,  albeit 
for  a  small  piece  of  the  macroscopic  system.  This  will  be  discussed  in  more 
detail  shortly. 

6. 1.2.  Slab  Models  These  models  are  the  preferred  geometry  in  electronic 
structure  methods  originating  from  solid  state  physics.  Periodic  boundary 
conditions  are  applied  to  a  unit  cell  representing  the  M/C  interface.  There¬ 
fore,  slab  models  are  restricted  to  coherent  interfaces,  which  means  that 
periodicity  parallel  to  the  interface  is  present  -  this  corresponds  to  a  “locked 
in”  interface  structure.  Of  course,  the  periodicity  may  have  a  long  repetition 
length,  which  may  not  be  modeled  properly  in  some  periodic  slab  calcula¬ 
tions.  There  are  two  types  of  slab  cells:  those  that  do  not  include  a  vacuum 
layer  and  those  that  do.  Dense  unit  cells  (Fig.3B)  physically  correspond  to 
sandwich  structures  with  two  M/C  interfaces  per  unit  cell.  From  an  ener¬ 
getic  point  of  view,  this  geometry  is  rather  restrictive,  because  it  requires 
sufficient  symmetry  such  that  the  interfaces  arc  identical  (in  order  for  the 
interface  energy  to  be  uniquely  defined).  Unfortunately,  this  requirement 
often  restricts  transverse  relaxation,  because  the  symmetry  “locks”  the  M/C 
interface.  On  the  other  hand,  if  symmetry  is  enforced,  this  unit  cell  has  no 
net  perpendicular  dipole  moment,  so  that  unphysical  electrostatic  coupling 
between  the  interfaces  is  avoided.  It  is  also  necessary  to  make  both  the 
metal  and  ceramic  layer  thick  enough  so  that  the  interfaces  do  not  interact 


16  Asbjorn  Christensen,  Emily  A.  Asche  ancl  Emily  A.  Carter 


via  electronic  structure  perturbations  or  strain,  in  order  to  have  such  a  model 
realistically  portray  a  single  M/C  interface. 

A  more  general  interface  geometry  is  shown  in  Fig.3C.  Physically  this 
corresponds  to  an  infinite  array  of  M/C  thin  film  couples  separated  by  vac¬ 
uum.  A  salient  point  is  that  the  vacuum  layers  should  be  thick  enough  that 
adjacent  M/C  slabs  do  not  interact.  Interaction  is  possible  in  two  ways: 
either  via  electronic  wavefunction  overlap  in  the  vacuum  or  via  Coulom- 
bic  multipoles.  The  former  interaction  is  usually  vanishing,  if  more  than 
~  10  A  of  vacuum  is  present.  The  latter  interaction  is  rather  long-ranged, 
but  fortunately  methods  have  been  devised  to  electrostatically  decouple  the 
slabs.108-110  Of  course,  it  is  required  that  both  the  metal  and  ceramic  layers 
arc  thick  enough  that  the  interface  and  surfaces  do  not  interact. 

The  slab  geometry  also  suffers  from  other  finite  size  effects.  If  the  extent 
of  the  unit  cell  parallel  to  the  interface  is  too  small,  artificial  strain  effects 
arc  introduced,  because  the  metal  and  ceramic  are  forced  to  be  coherent 
by  the  periodic  boundary  conditions.  Of  course,  this  may  be  eliminated 
by  enlarging  the  unit  cell,  which  unfortunately  leads  to  very  computer¬ 
intensive  calculations,  as  is  the  case  with  the  cluster  models.  However  for 
the  slab  model,  the  oscillations  in  the  electronic  density  of  states  arc  not  as 
dramatic  when  varying  the  number  of  atoms  as  in  the  case  with  clusters. 
This  is  because  the  slab  is  infinite  parallel  to  the  interface.  This  implies  the 
spectrum  is  continuous,  and  the  metal  slab  does  not  have  an  artificial  band 
gap,  unlike  the  metal  cluster. 

The  artificial  requirement  of  the  unit  cell  periodicity  perpendicular  to 
the  interface  may  be  avoided,  if  a  Green’s  function  technique  is  used.  Such 
techniques  have  been  used  for  metal-metal  interfaces1 1 1  with  high  symmetry, 
but  no  such  calculations  for  M/C  interfaces  have  been  reported  yet  to  our 
knowledge. 

6.1.3.  Interface  Stoichiometry  If  no  experimental  evidence  is  available 
concerning  the  orientation  of  the  metal  and  ceramic  crystals  with  respect 
to  each  other  for  a  coherent  interface,  one  needs  somehow  to  produce  a 
reasonable  guess  as  to  how  the  parent  crystals  may  match  up  at  the  interface. 
The  possibilities  arc  immense,  especially  if  the  ceramic  has  low  symmetry: 
first  one  needs  to  consider  which  faces  from  each  parent  crystal  will  match 
up.  Symmetry  considerations  arc  often  helpful  at  this  point.  Free  energy 
consideration  might  also  be  guiding:  weakly  interacting  M/C  pairs  arc  likely 
to  match  up  on  their  low  energy  surfaces.  Strongly  interacting  M/C  pairs 
arc  more  likely  to  match  up  on  their  high  energy  surfaces,  because  these 
generally  contain  more  unsaturated  atoms  ready  for  bonding.  However,  if 
a  new  reaction  phase  forms  at  the  interface,  the  situation  becomes  more 
complicated.  If  the  bonding  mechanism  is  dominantly  electrostatic,  image 
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theory,  which  we  return  to  later,  predicts  the  most  polar  ceramic  surface  to 
bind  most  tightly  to  the  metal,  i.e.,  it  has  the  largest  work  of  separation  W. 
However,  this  does  not  mean  that  polar  interfaces  arc  most  stable  from  a 
thermodynamical  point  of  view;  the  most  stable  M/C  interface  that  can  be 
formed  between  the  pair  M  and  C  is  the  one  with  smallest  free  energy,  i.e. 
smallest  interface  tension  yMC: 

yMC  =  yM  +  yC-  W~Y_  RANi  (1) 

i 

This  is  the  Dupre  relation,  generalized  to  allow  mass  exchange,  e.g.,  with  the 
ambient  gas  or  with  the  bulk,  at  chemical  potentials  jy.  AN  ^  is  the  change  in 
abundance  of  species  i  in  the  interface  region  on  forming  the  interface,  and 
yM,YC  arc  surface  tensions  of  hare  metal  and  ceramic,  respectively.  Thus, 
because  a  polar  surface  is  energetically  very  unfavorable  (theoretically  yc 
diverges  for  a  perfect,  infinite  polar  ceramic  surface),  it  may  not  be  favorable 
to  form  an  interface  with  a  polar  ceramic  termination,  even  though  the 
interface  binding  is  relatively  strong.  Another  approach  for  guessing  the 
ceramic  termination  at  the  interface  is  to  estimate  anion/cation/metal  bond 
combinations  from  the  appropriate  bulk  phases.  But  no  foolproof  rules 
can  be  formulated  on  how  the  ceramic  is  terminated  at  the  M/C  interface. 
Stoichiometry  is  a  complicating  aspect  in  this  context,  because  it  depends 
on  the  chemical  potentials  present;  and  the  interface  chemistry  changes 
dramatically  with  stoichiometry,  as  we  will  discuss  later. 

6.1.4.  Lattice  Misfit  Having  settled  on  some  metal  and  ceramic  surfaces 
that  we  think  will  match,  we  must  determine  the  relative  orientation  and 
translation  of  these  surfaces  with  respect  to  each  other.  Certain  directions  in 
the  metal  and  ceramic  arc  likely  to  be  aligned  for  a  stable  interface.  There  will 
often  be  multiple  minima,  corresponding  to  different  lock-in  possibilities  for 
the  coherent  interface.  Lastly,  the  size  and  shape  of  the  interface  unit  cell 
needs  to  be  determined,  if  we  assume  a  coherent  interface,  which  is  implicit 
if  periodic  boundary  conditions  arc  applied.  A  realistic  unit  cell  will  of 
course  correspond  to  low  strain  on  both  the  metal  and  ceramic  side. 

This  vast  number  of  possibilities  calls  for  a  systematic  procedure  to  iden¬ 
tify  a  subset  of  the  most  likely  interface  matchings  of  the  parent  crystals.  This 
subset  will  then  be  the  starting  point  for  atomistic  modeling.  The  question 
about  unit  cell  size  and  shape  is  relatively  simple  to  address.  Many  related 
procedures  based  on  linear  elasticity  theory  and  lattice  strain  estimates  may 
be  adopted.  The  basic  situation  is  sketched  in  Fig.  4:  an  overlayer  unit  cell 
A  needs  to  be  matched  together  with  a  substrate  unit  cell  B.  Matching  pairs 
of  unit  cells  are,  in  general,  multiples  of  primitive  cells  in  the  interface  plane 
for  the  metal  and  ceramic,  respectively. 
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Figure  4  Schematic  picture  of  M/C  interface  lattice  matching. 

A  :  Arbitrary  overlayer  unit  cell  constituted  of  many  primitive  cells  spanned  by  the  vectors 
cti ,  &2  • 

B  :  Arbitrary  substrate  unit  cell  constituted  of  many  primitive  cells  spanned  by  the  vectors 
bi  ,b2. 

T  :  Linear  transformation  between  A  and  B. 

A+B  :  The  “overlap”  Cl  between  cells  A  and  B. 


The  most  crude  approach  to  quantify  the  commensurability  of  the  two 
given  unit  cells  is  to  assign  a  mismatch  factor  based  on  the  overlap  between 
the  cells  A  and  B,  i.e. 


q  =  l 


m\ 

IAI  +  IBI 


(2) 


where  jA|,  |B|  and  jOj  signify  the  area  of  cell  A,  B  and  the  overlap  area, 
respectively.  This  quantity  vanishes,  if  A  =  B  and  increases  with  decreasing 
overlap  Cl,  as  sketched  in  Fig.  4.  Thus  one  would  like  to  choose  an  interface 
matching  that  minimizes  p. 

Another  geometrical  approach  was  taken  by  Bolding  and  Carter,112  who 
considered  the  matching  of  orthogonal  unit  cells.  By  taking  multiples  of 
each  unit  cell  a  suitable  number  of  times  in  each  direction,  an  arbitrarily 
low  misfit  might  be  achieved,  in  the  same  way  that  an  arbitrary  irrational 
number  might  be  approximated  by  a  rational  number,  if  the  numerator  and 
denominator  are  sufficiently  large.  More  specifically,  they  defined  the  strain 
variables 


6 


nl 

1 


lai  —  Tibi 
ubi 


(3) 


Atomic-level  properties  of  thermal  barrier  coatings :  characterization  of  metal-ceramic  interf races  19 


k  _  ka2  -  mb2 
2  mb2 

where  ai ,  a2,  hi ,  b2  arc  the  length  of  the  basis  vectors  spanning  the  primitive 
unit  cells  of  which  A  and  B  are  multiples,  see  Fig.  4.  (n,  m,  l,  k)  are  integers 
giving  the  unit  cell  enlargement  along  each  basis  vector.  Then  Bolding  and 
Carter  proposed  minimizing  the  quantity 

^?llc=w1|5?l|+w2|65llcl  (4) 

which  obviously  is  positive  definite  and  vanishes  at  perfect  coherency  be¬ 
tween  A  and  B.  Here  Wi ,  w2  arc  weight  factors  that  might  be  related  to  the 
elastic  properties  of  metal  and  ceramic.  In  Bolding  and  Carter’s  case,  these 
weight  factors  were  set  to  unity. 

A  more  general  approach  would  be  to  directly  consider  the  elastic  energy 
associated  by  deforming  cell  A  into  cell  B  so  as  to  achieve  coherency.  If  this 
deformation  is  designated  by  the  2x2  linear  transformation  matrix  T 

B  =  TA  (5) 

where  the  matrices  A  and  B  contain  the  vectors  spanning  the  (almost)  match¬ 
ing  pair  of  metal  and  ceramic  unit  cells,  it  is  easy  to  show  that  the  corre¬ 
sponding  exact  strain  tensor  e  and  the  corresponding  elastic  energy  u  is 
given  by 

e  =  ^[T'T-I]  (6) 

1 

u  =  -e:y:e  (7) 

It  is  important  to  use  the  exact  strain  tensor  definition,  Eq.  (6),  to  achieve 
rotational  invariance  with  respect  to  lattice  rotation;  the  conventional  linear 
strain  tensor  only  provides  differential  rotational  invariance  of  u  in  Eq. 
(7)  ii3, 114  a  hierarchy  of  approximations  may  be  used  for  the  elastic  tensor 
y.  The  most  rigorous  approach  is  to  transform  the  bulk  elastic  tensor  c 
according  to 

Ypqrs  =  Rpi^qjCi.j]cxRlcrkls  (8) 

where  a  summation  convention  over  same  indices  applies  and  R  is  the  rotation 
matrix  from  interface  to  bulk  crystal  Cartesian  frames.113  However,  it  is  not 
guaranteed  that  a  thin  film  has  the  same  ratio  between  elastic  constants  as 
the  bulk  crystal.  Generally,  six  elastic  constants  are  inequivalent  in  2D 
elasticity  theory,  which  follows  from  the  permutational  symmetry  of  c.  As 
an  approximation  we  may  assume  only  Cn,  Cu  and  C44  arc  nonvanishing, 
which  is  true  for  cubic  symmetry.  The  lowest  level  of  approximation  is  to 
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assume  isotropy,  in  which  case  the  relation  Cu  =  C12  +  2C44  applies  and 
only  a  ratio,  like  C12/C1 1 ,  needs  to  be  estimated. 

The  approach  assumes  a  fixed  substrate;  the  most  sophisticated  procedure 
would  be  to  minimize  the  elastic  energy  u  =  usubstTQte  +  uoverlaver_  This 
properly  results  in  a  small  overall  rescaling  of  the  elastic  energies.  However, 
due  to  the  physical  simplicity  of  this  elastic  model,  it  is  doubtful  whether  a 
better  description  is  obtained  unless  both  substrate  and  overlayer  arc  highly 
anisotropic.  Further,  if  material  A  is  grown  onto  a  bulk  substrate  B,  it  is 
unlikely  that  a  transversal  deformation  of  the  substrate  would  take  place  in 
reality  anyway. 

It  is  instructive  to  compare  this  model  to  that  of  Bolding  and  Carter112 
discussed  above,  for  the  case  where  both  substrate  and  overlayer  have  or¬ 
thogonal  lattices  (in  the  interface  plane).  If  we  multiply  the  overlayer  unit 
cell  (l,lc)  times  along  cij ,  0.2  and  correspondingly  the  substrate  unit  cell 
(u,  m)  times  along  bi ,  b2,  the  transformation  matrix  T  in  Eq.(5)  becomes 

/  nbj_  a  \ 

ynmlk  =  ^  U,  ^  \  (9) 

Inserting  this  into  the  strain  tensor  Eq.(6),  we  get  the  harmonic  energy, 
Eq.(7),  to  second  order  in  the  strain  variables  of  Bolding  and  Carter, 

Uumtic  =  ^  +  1T22(6^)2  +  y126?l62*  (10) 

Compared  to  the  misfit  factor  of  Bolding  and  Carter  uj(jTlk,  Eq.(4),  this 
expression  contains  a  cross  term  and  scales  differently  in  the  strain  variables. 
Both  expressions  will  be  zero  for  perfect  coherency  (61  =  62  =  0),  but  for 
competing  interface  lock-in  possibilities,  one  might  expect  Eq.(10)  to  be 
qualitatively  more  accurate.  In  conclusion,  to  find  the  most  favorable  unit 
cell  lock-ins  for  a  given  substrate  and  overlayer  lattice,  from  an  elasticity 
point  of  view,  we  simply  need  to  scan  an  expression  like  Eq.(10)  in  four 
indices  (u,  m,  l,  k)  to  find  a  subset  of  likely  coherent  interface  unit  cells. 
This  is  easily  done  on  a  computer. 

In  addition  to  the  elastic  energy,  of  course,  the  chemical  interaction  en¬ 
ergy  between  metal  and  ceramic  must  be  accounted  for.  The  possibility 
of  competition  between  elasticity  and  chemistry  exists  such  that  a  rather 
elastically  strained  interface  combination  may  be  more  stable  than  an  un¬ 
strained  one,  because  the  strained  one  may  have  favorable  chemical  bonding 
between  metal  and  ceramic  atoms  at  the  interface.  Therefore,  other  interface 
combinations  than  the  least  strained  need  to  be  considered  in  general. 

We  now  move  on  to  discuss  how  one  can  assess  structure  and  energetics 
of  such  interfaces  from  theory,  starting  from  the  most  accurate  models  and 
ending  with  the  most  approximate. 
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6.2.  Ab  Initio  Techniques 

First  principles  techniques  can  be  distinguished  from  semiempirical  meth¬ 
ods  in  that  they  include  quantum  mechanical  effects  explicitly  using  exact 
or  moderate  approximations  for  electronic  exchange  and  correlation.  Gen¬ 
erally,  explicit  inclusion  of  quantum  mechanical  effects  is  computationally 
demanding,  which  restricts  ab  initio  methods  to  small  systems,  typically 
less  than  50  inequivalent  atoms  on  modern  workstations  (1999).  Using  mas¬ 
sively  parallel  computers,  up  to  300  inequivalent  main  group  atoms  can  be 
handled  ab  initio  today.*  In  specialized  cases,  significantly  more  atoms  can 
be  handled  by  so-called  O  (N )  ab  initio  techniques,  where  the  computational 
load  scales  essentially  linearly  with  the  number  of  atoms  N.  For  insulating 
systems,  tight  binding  calculations  have  shown  that  electronic  density  matrix 
methods1 15  112  may  be  applied  since  the  electronic  density  matrix  has  a  finite 
range,  and  simulations  with  650  atoms  on  a  modern  workstation  have  been 
reported. 119  Alternatively,  divide  and  conquer120  and  localized  wave  function 
approaches121-124  arc  well  suited  to  achieve  O  (N )  scaling  for  insulating  sys¬ 
tems.  If  an  artificial  nonzero  temperature  is  assumed  for  the  electron  system, 
density  matrix117, 125  127  and  localized  wave  function  approaches128  arc  also 
feasible  for  metallic  systems,  but  imposing  a  nonphysical  electronic  temper¬ 
ature  -  unless  it  is  removed  by  the  end  of  the  calculation127  -  is  a  dubious 
approximation  for  systems  where  the  electronic  density  of  states  has  much 
structure  around  the  Fermi  level,  as  in  the  case  for  transition  metal/ceramic 
interfaces.  For  free -electron-like  bulk  systems.  Orbital-free  Density  Func¬ 
tional  Theory  involving  kinetic  energy  density  functionals  has  allowed  the 
treatment  of  ~6000  symmetry-inequivalent  metal  atoms.129  The  drawback 
of  this  method  is  that  it  is  limited  by  the  accuracy  of  the  kinetic  energy 
density  functional,  which  is  not  known  exactly  for  many-electron  systems. 
Improving  such  functionals  is,  however,  an  active  area  of  research.130-144 

The  restriction  on  the  number  of  inequivalent  atoms  implies  that  only 
properties  involving  short  length  scales  can  be  treated  accurately,  unless 
these  properties  arc  consistent  with  periodic  boundary  conditions.  This 
restriction  is  rather  severe  in  the  case  of  a  M/C  interface,  because  the  real 
interface  is  often  transversely  aperiodic  and  may  contain  long-ranged  strain 
fields  and  defect  structures  perpendicular  to  the  interface,  which  will  require 
an  extremely  large  unit  cell  to  model.  Thus  only  idealized  M/C  interfaces 
having  small  unit  cells  may  be  treated  realistically  by  ab  initio  methods. 


t  The  definition  of  the  concept  ab  initio  is  somewhat  fluid,  but  usually  means  essentially  free  of  empirical 
parameters.  Being  ab  initio  is  generally  considered  a  quality  stamp,  but  one  should  be  aware  that  some 
approaches,  which  are  formally  parameter-free,  e.g.  Thomas-Fermi  density  functional  theory  or  ab  initio 
tight  binding,  fail  miserably  outside  limited  application  ranges. 
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This  excludes  many  M/C  combinations  with  large  lattice  mismatch  between 
the  metal  and  the  ceramic,  because  a  small  interface  unit  cell  implies  that 
either  the  metal  or  the  ceramic  must  be  strained  to  an  unphysical  extent, 
leading  to  unphysical  results  induced  by  the  artificial  periodic  boundary 
conditions.  Still,  for  M/C  interfaces  where  the  natural  lattice  mismatch 
is  low,  ab  initio  studies  of  idealized  M/C  interfaces  arc  instructive,  giving 
valuable  insights  into  quantum  effects  at  the  interface.  We  now  discuss  what 
has  been  accomplished  with  ab  initio  calculations  of  idealized  interfaces. 

6.2.1.  Quantum  Chemical  Approaches  A  more  precise  description  for  this 
class  is  full  wavefunction  methods,  where  the  basic  variable  is  the  full  many- 
body  wavefunction.  The  main  problem  with  full  wavefunction  approaches  is 
that  the  computational  load  increases  drastically  with  the  number  of  electrons 
N.  At  the  Hartree-Fock  level,  the  load  increases  as  N3-4,  and  the  scaling 
with  N  increases  steadily,  the  more  complete  the  inclusion  of  electronic 
correlation.  Since  realistic  modeling  of  interfaces  requires  quite  large  unit 
cells  or  large  clusters,  very  few  full  wavefunction  studies  at  and  beyond  the 
Hartree-Fock  level  have  been  reported  in  the  literature.  The  interaction  of 
a  Cu145  and  a  Pd146  atom  with  a  single  MgO  molecule  have  been  studied 
using  multireference  configuration  interaction.  In  both  cases,  binding  of  the 
metal  to  the  O-atom  in  the  MgO  molecule  was  predicted  as  most  stable  (i.e. 
a  lineal-  molecule),  with  an  “adsorption”  energy  of  2.28  eV  and  2.65  eV  for 
Cu  and  Pd,  respectively.  This  is  in  accordance  with  the  trend,  observed  by 
less  accurate  electronic  structure  methods,  that  metal  atoms  adsorb  on  top 
of  the  O  atoms  on  the  MgO(OOl)  surface. 

Hartree-Fock  studies  have  been  undertaken  for  cluster  models  of  the 
Cu/MgO147  and  Ni/Al203148  interfaces.  Both  studies  find  significant  charge 
donation  from  the  metal  atom  to  the  oxide  conduction  states.  Hartree-Fock 
with  a  posteori  electron  correlation  taken  from  a  free  electron  gas  (HF- 
CC)154  in  the  supercell  approach  has  been  applied  to  Ag/a-Al203155  and 
Ag/MgO(001).155- 156  Ag  was  found  to  be  weakly  bound  to  the  completely 
O-terminated  a-Al203(0001)  surface,  although  significant  charge  transfer 
occurred.  In  agreement  with  density  functional  calculations,  HF-CC  calcu¬ 
lations  find  Ag  physisorbed  above  the  O  atoms  in  the  MgO(OOl)  surface. 
The  adsorption  energy  above  the  O(Mg)  site  is  0.47(0.11)  J/m2,  in  good 
agreement  with  the  experimental  value,  0.45  J/m2  for  a  thick  Ag  film.157 
However,  this  number  is  likely  to  be  an  average  of  adsorption  of  Ag  over 
both  O  and  Mg  in  separate  domains,  due  to  the  (small)  lattice  misfit  be¬ 
tween  Ag  and  MgO.  A  principal  problem  with  the  HF-CC  method  is  that 
the  a  posteori  free  electron  gas  electronic  correlation  is  not  consistent  with 
the  exact  exchange  energy  used  in  the  underlying  Hartree-Fock  method  - 
important  error  cancellations  between  exchange  and  correlation  are  lost.158 
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Initial  tests154  indicated  a  confidence  level  comparable  to  the  best  density 
functional  models  for  exchange  and  correlation,  but  the  approach  needs  to  be 
tested  on  more  systems  to  see  if  systematic  errors  exist.  Related  in  spirit  to 
the  HF-CC  method  is  the  B3LYP159  approach,  which  evaluates  the  electronic 
exchange-correlation  energy  as  a  weighted  average  of  the  exact  exchange 
energy  and  the  electron  gas  results  for  the  exchange-correlation  energy.  The 
weighting  parameters  for  this  average  were  determined  semiempirically, 
which  makes  this  approach  aesthetically  unpleasant.  However,  for  many 
molecules  and  finite  clusters  this  approach  has  delivered  accurate  results. 
An  instructive  study  has  appeared,160  which  compares  the  adhesion  of  Cu 
atoms  on  MgO(OOl)  using  many  electronic  structure  methods,  ranging  from 
various  forms  of  density  functional  theory  to  post-Hartree-Fock  approaches. 
Here,  the  results  of  the  B3LYP  method  agreed  rather  well  with  the  results 
obtained  by  the  best  post-Hartree-Fock  approach,  which  estimated  the  bind¬ 
ing  energy  of  the  Cu  atom  on  MgO(OOl)  to  0.40  ±  0.05  eV.  This  is  probably 
the  most  accurate  theoretical  estimate  of  this  quantity  today.  This  paper 
also  illustrated  quite  dramatically  the  sensitivity  of  the  results  to  choice  of 
DFT  functional,  which  should  sound  a  note  of  caution  regarding  the  general 
applicability  and  transferability  of  current  implementations  of  DFT. 

In  summary,  full  wavefunction  methods  are  the  most  accurate  ab  initio 
methods,  if  electronic  correlation  is  accounted  for  properly  and  the  basis 
set  for  expansion  of  the  electronic  wavefunction  is  sufficiently  complete  - 
this  might  be  difficult  to  judge  by  a  non-expert.  Furthermore,  if  a  cluster 
model  is  used,  it  needs  to  be  of  sufficient  size  or  be  properly  embedded  for 
conclusions  to  be  representative  of  an  extended  M/C  interface. 

6.2.2.  Density  Functional  Theory  For  condensed  phase  systems.  Density 
Functional  Theory  (DFT)158, 161- 162  methods  constitute  the  optimal  compro¬ 
mise  between  accuracy  and  efficiency  of  all  ab  initio  methods  available 
today.  The  key  point  of  DFT  is  to  show  that  the  exact  quantum  mechan¬ 
ical  total  energy  E  is  a  functional  only  of  the  total  electron  density  p(r), 
a  3-dimensional  function,  which  corresponds  to  the  intractable  many-body 
electronic  wavefunction  W,  a  3/V-dimcnsional  function  where  N  is  the  num¬ 
ber  of  electrons.  The  total  energy  within  DFT,  E[p],  is  related  to  the  total 
energy  expressed  in  terms  of  T,  E(T),  as 

E[¥]  =  <  W\telec  +  'QjV  >  +Tion 

=  Tselec[p]  +  Uc|ass-ica;[p]  +  uxc[p]  +  Tlon  (11) 
=  E[p] 

In  principle,  DFT  is  thus  a  formally  exact  mean  field  theory.  Here  Tlon  is  the 
ionic  kinetic  energy,  UC|QSS|cax  is  the  classical  mean  field  Coulomb  energy 
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for  ionic  cores  and  charge  distribution  p,  T®lec  represents  the  kinetic  energy 
of  the  (noninteracting)  electrons,  and  Uxc  is  a  term  describing  lowering  of 
electron-electron  repulsion  by  correlated  electronic  motion  and  electronic 
exchange,  as  well  as  a  kinetic  energy  component  due  to  electron-electron 
interactions.  Uxc  in  Eq.  (1 1)  is  the  only  term  whose  functional  form  is  not 
known  exactly:  the  most  commonly  used  approximations  to  it  today  arc  the 
Local  Density  Approximation  (LDA)  and  the  Generalized  Gradient  Approx¬ 
imation!  s)  (GGA).149  A  variety  of  slightly  different  GGA  parameterizations 
have  been  proposed  over  the  years;150-153  indeed  the  choice  of  which  GGA 
functional  to  use  is  a  nonsystematic  aspect  of  DFT  calculations. 

A  slightly  lower  level  of  approximation  than  the  LDA  is  the  Xa  method, 
where  the  correlation  energy  is  assumed  to  be  proportional  to  the  exchange 
energy,  for  which  the  LDA  uniform  electron  gas  expression  is  used.  The 
Xa-approximation  has  been  used  for  many  pioneering  metal/ceramic  model 
studies.  A  consensus  has  developed  in  the  literature154, 160,163-165  that  the 
GGA  on  average  is  better  than  the  LDA.  But  in  certain  cases,  the  GGA  tends 
to  overcorrect  the  LDA. 166, 167  For  ionic  surfaces,  a  significant  lowering  of  the 
surface  energy  going  from  LDA  to  GGA  has  been  noticed  in  some  cases.168 
Another  major  drawback  of  the  LDA/GGA  approximations  in  this  context  is 
that  they  are  not  able  to  describe  strongly  correlated  oxides,  such  as  many  3d 
transition  metal  oxides.169  These  cases  with  breakdown  of  the  LDA/GGA 
arc  due  to  two  shortcomings  of  the  LDA/GGA.  The  first  is  the  general 
underestimation  of  the  band  gap,  which  can  be  traced  back  to  a  discontinuity 
feature  in  the  exchange-correlation  potential  as  function  of  electron  filling162 
not  captured  by  the  LDA/GGA.  The  second  is  the  self-interaction  problem 
in  the  LDA/GGA  in  which  electrons  artificially  interact  with  themselves  due 
to  the  electron  gas  approximation  for  the  exchange-correlation  energy.  This 
may  be  remedied  by  considering  a  self-interaction  correction  (SIC)170, 171 
scheme  to  the  LDA/GGA,  e.g.  the  LDA+U  method.172 

For  insulating  materials,  such  as  ceramics,  this  normally  does  not  pose 
problems,  because  the  conduction  band  is  empty  and  the  band  gap  error  is 
not  reflected  energetically.  However,  for  M/C  interfaces  there  may  be  some 
reason  to  worry,  since  mixing  between  oxide  conduction  bands  and  metal 
states  will  influence  the  magnitude  of  adhesion  and  charge  transfer.  These 
states  may  become  filled,  depending  on  the  Fermi  level  on  the  metal  side. 
Further  research  on  this  issue  is  warranted. 

Because  DFT-based  techniques  have  the  electronic  density  p(r)  as  the 
basic  variable,  the  computational  load  scales  moderately  with  the  number  of 
electrons  IN ,  O  ( N 1  _3 ) .  Thus  they  are  able  to  handle  significantly  more  atoms 
than  traditional  quantum  chemical  approaches  that  retain  the  full  electronic 
many-body  wavefunction  as  the  basic  variational  quantity.  This  favorable 


Atomic-level  properties  of  thermal  barrier  coatings :  characterization  of  metal-ceramic  interfraces  25 


scaling  currently  makes  DFT-based  techniques  most  promising  for  ab  initio 
studies  of  M/C  interfaces,  and  therefore  we  will  emphasize  this  group  of 
methods  in  our  review. 

6.2.2. 1.  Applications  of  D FT  to  M/C  Interfaces  An  increasing  number 
of  studies  using  self-consistent  DFT-based  techniques  have  been  reported 
recently.  The  most  popular  M/C  interface  among  theoreticians  has  been 
Ag/MgO(001),  due  to  its  small  unit  cell  and  small  lattice  mismatch  of 
3%  and  hence  simple  epitaxial  character.  Furthermore,  this  interface  is 
well  characterized  experimentally.  Generally,  studies  based  on  electronic 
structure  methods  agree  that  this  interface  is  characterized  by  weak  phys¬ 
ical  rather  than  strong  chemical  adhesion,  with  marginal  charge  transfer 
to  Ag,  and  the  O-site  as  most  favorable  adsorption  site.155, 156,173-175  The 
MgO(OOl)  surface  is  stoichiometric  (i.e.,  there  is  an  equal  number  of  Mg 
and  O  ions  on  the  surface)  and  O  is  nominally  in  the  02~  charge  state.  As 
also  noted  by  Finnis,94  it  is  interesting  to  observe  the  scattering  in  adhe¬ 
sion  energies  obtained  by  different  basis  set  expansions.  For  Ag  adsorbed 
on  top  of  O(Mg),  the  Full-Potential  Linear  Muffin  Tin  Orbital  (FP-LMTO) 
method  gives  1.59(0.78)  J/m2  for  3  layers  of  Ag,176  the  Full-Potential  Linear 
Augmented  Plane  Wave  (FLAPW)  method  gives  0.53(0.53)  J/m2  for  1  Ag 
layer,174  whereas  a  Gaussian-based  Linear  Combination  of  Atomic  Orbitals 
representation  gives  1.9(1.08)  J/m2  for  1  Ag  layer.177  All  of  these  studies 
were  performed  within  the  LDA,  which  is  notorious  for  its  tendency  towards 
overbinding.158  This  overbinding  trend  is  followed  for  the  Ag/MgO  system, 
confirmed  by  comparing  to  the  experimental  value  0.45  J/m2  for  a  thick  Ag 
film.157  Again,  this  value  is  likely  to  be  an  average  of  alternating  domains 
with  Ag  adsorbed  over  O  and  Mg  respectively,  due  to  misfit  dislocations  in¬ 
duced  by  strain.  The  scatter  in  the  theoretical  results  is  not  explained  by  the 
extent  to  which  relaxations  in  the  interface  region  were  included;  all  authors 
relaxed  the  Ag-0  distance,  and  none  included  full  relaxations  in  the  oxide. 
The  FP-LMTO  study176  relaxed  the  oxide  planes  perpendicular  to  the  sur¬ 
face,  the  FLAPW  study  kept  the  oxide  frozen.  The  free  MgO(OOl)  surface 
termination  displays  a  rumpling  of  the  ions;  however,  our  own  DFT  calcula¬ 
tions,178-  179  which  feature  full  relaxation,  show  that  the  Ag-film  smooths  the 
MgO(OOl)  surface  rumpling,  so  that  the  oxide  layer  nearest  to  the  interface 
is  flat.  Furthermore,  the  oxide  interlayer  distances  are  as  in  bulk  MgO  to 
within  0.005  A.  Fig.  5  displays  a  repeated  unit  cell  from  our  calculations 
of  an  Ag  layer  on  top  of  a  MgO(OOl)  surface.  Thus,  in  this  particular  case, 
the  structural  constraint  in  the  FLAPW  study  did  not  affect  the  results.  The 
large  variations  in  theoretical  predictions  remains  an  unsatisfying  aspect  of 
these  calculations. 
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Figure  5  Ag/Mg0(001)  charge  density  plot  cross-sectional  view. 

There  are  seven  Ag  sites  on  top  of  the  layer  of  seven  O  sites.  The  Mg,  with  little  valence 
charge,  appear  faintly  beneath  the  O  sites.  Another  O  layer  and  Mg  layer  appears  at  the 
bottom.  The  contour  lines  at  the  top  of  the  box  are  a  result  of  periodic  boundary  conditions. 


Polar  M/C  interfaces  have  only  been  studied  sparsely.  The  polar  (111) 
and  nonpolar  (100)  Cu/MgO  interfaces  were  compared  in  a  slab  calculation 
within  the  LDA. 180  As  expected  from  the  image  model  theory,  the  polar 
interface  displayed  a  considerably  higher  work  of  separation.  In  addition 
these  authors  found  a  larger  charge  transfer  and  Cu-MgO  orbital  mixing  for 
the  polar  interface. 
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Generally,  ah  initio  methods  have  found  that  metals  adsorb  at  the  O-site 
of  MgO(OOl)  and  that  the  Mg-site  corresponds  to  a  local  maximum  on  the 
adsorption  potential  energy  surface  (PES),  but  the  accord  ends  here.  For 
Ag/MgO(001),  Hong  et  al . 177  and  Schonberger  et  al.1 76  find  a  corrugation 
of  the  PES  of  the  order  of  50  %  of  the  cohesive  energy  whereas  Li  et  al. 174 
find  an  extremely  flat  PES  with  corrugation  less  than  1  %.  For  comparison, 
for  Pd/MgO(001),  Goniakowski181  finds  a  corrugation  of  order  20  %.  Thus 
even  though  these  systems  appear-  easy  to  study,  the  disagreement  upon  even 
qualitative  features  of  the  interaction  is  rather  disturbing. 

It  is  an  interesting  question  whether  the  GGA  improves  the  agreement  with 
experiments,  as  is  often  the  case,  or  even  changes  trends  predicted  within  the 
LDA.  We  investigated  this  by  calculating  the  Ag/MgO(001)  adhesion  within 
both  the  LDA  and  GGA,  using  state-of-the-art  ultrasoft  pseudopotentials178 
in  a  converged  planewave  basis.  For  Ag  adsorbed  on  top  of  O(Mg),  we 
obtain  adhesion  energies  of  0.62(0.27)  J/m2  within  the  LDA  and  0.23(0.06) 
J/m2  within  the  GGA.  This  indicates  that  the  GGA  overcorrects  the  LDA 
somewhat,  compared  to  the  experimental  value  of  0.45  J/m2.157  Also,  a 
sensitivity  on  details  of  how  the  GGA  functional  is  parameterized  has  been 
noticed.160-182 

Only  recently  have  other  M/C  interface  studies  appeared  that  tested 
the  GGA  in  this  context.  Pacchioni  and  Rosch183  studied  Ni  and  Cu  on 
MgO(OOl)  in  a  cluster  model.  As  might  be  expected,  they  found  that  Ni 
binds  more  strongly  than  Cu,  due  to  the  open  d-shell  of  Ni.  In  both  cases, 
the  O-site  is  favored.  They  found  a  significant  lowering  of  the  adsorption 
energy  for  Cu  when  using  the  GGA,  0.65  J/m2,  compared  to  the  value  2.54 
J/m2  obtained  recently  in  a  cluster  calculation  using  the  LDA.184  Pacchioni 
and  Rosch  also  noted  that  in  these  cases  metal-metal  bonds  were  stronger 
than  metal-substrate  bonds,  thus  predicting  a  3D  (cluster)  growth  mode  as 
opposed  to  layer-by-layer  growth  (wetting).  A  more  systematic  study  of 
transition  metal  adsorption  on  MgO(OOl)  using  GGA  was  presented  by  Yu- 
danov  et  al.}15  using  a  cluster  model.  Interestingly,  they  found  no  obvious 
correlation  between  cohesion  and  d-band  filling  as  indicated  by  other  stud¬ 
ies.  They  found  that  single  Cu,  Ag,  Au,  Cr,  and  Mo  atoms  exhibited  weak 
interface  bonds,  whereas  Ni,  Pd,  Pt  and  W  atoms  formed  strong  bonds.  For 
Ag/MgO,  they  obtained  0.36  J/m2,  similar  to  our  findings  for  a  coherent  Ag 
monolayer. 

Square -planar  M4  cluster  adsorption  for  M  =  (Ni,  Cu,  Pd,  Ag)onMgO(OOl) 
has  been  studied  within  the  GGA. 1 85  Although  here  open  shell  transition  met¬ 
als  adsorbed  stronger,  metal-metal  and  metal-oxide  bond  competition  were 
found  to  be  complex.  An  extensive  study  of  Cu  growth  on  MgO(OOl)  was 
performed  for  clusters  with  1-13  Cu  atoms.186  Here  again,  the  authors  found 
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the  metal-metal  bonds  stronger  than  the  metal-ceramic  bonds,  thus  favoring 
a  3D  (cluster)  growth  mode.  Relatively  weak  cohesion  for  the  Cu/MgO 
interface  is  also  expected  on  the  basis  of  the  large  lattice  mismatch  (13  %) 
between  the  metal  and  ceramic. 

Electronic  structure  calculations  of  metal/Al203  interfaces  were  initiated 
by  Johnson  and  Pepper.187  They  modeled  this  interface  with  a  metal-AlOg9- 
cluster  using  the  Xa-method,  with  the  metal  atom  adsorbed  in  an  O-triangle 
emulating  a  hollow  site  on  a  02~ -terminated  surface.  This  study  concluded 
the  bonding  became  weaker  as  the  transition-metal  adsorbate  became  more 
noble,  which  was  explained  in  the  conventional  orbital  mixing  picture,  where 
antibonding  levels  become  filled  in  the  right  side  of  the  transition  series.  This 
picture  was  essentially  confirmed  by  later  investigators.  Most  authors  focus 
on  the  O-terminated  <x— A^OafOOOl)  surface,  which  is  nonstoichiometric,  as 
the  substrate  for  metallic  growth.  Under  oxygen  deficient  conditions,  how¬ 
ever,  it  is  not  clear  that  the  O-termination  is  appropriate  for  growth  modeling. 
Grazing  incidence  X-ray  scattering  data  indicate  that  the  <x— Al203(0001)  is 
Al-terminated  under  ultrahigh  vacuum  (UHV).188- 189  Furthermore,  DFT  cal¬ 
culations  predict  that  the  charge-neutral  O-terminated  <x— Al203(0001)  has  a 
significantly  higher  cleavage  energy190  than  Al-terminated  <x— AUOsfOOOl). 
Experimentally,  when  Ni  was  deposited  on  <x— AI2O3  in  UHV,  no  Ni-0  inter¬ 
action  was  detected  by  XPS.191  It  is  certainly  reasonable  to  consider  Ni-Al 
bonds,  since  NiAl  is  a  very  stable  alloy.  In  any  case,  discussions  about 
which  Al203-substrate  termination  at  the  clean  metal/ AI2O3  interface  is 
most  relevant  can  be  somewhat  academic,  since  reaction  phases,  like  spinels 
(MAI2O4),192  arc  often  formed  at  the  interface.  Our  point  of  view  is  that 
most  structural  possibilities,  at  this  stage,  are  interesting,  since  they  help 
to  furnish  a  general  understanding  of  the  electronic  structure  of  the  M/C 
interface  under  various  conditions. 

We  have  used  DFT  to  study  the  N1/AI2O3  interface.  These  calcula¬ 
tions  were  performed  within  the  planewave  DFT  code  VASP  using  ultra- 
soft  pseudopotentials.178’ 193  Calculations  with  one  and  three  layers  of  <x- 
AI2O3  deposited  onto  the  Ni(lll)  surface  were  performed.  Both  A1  and 
O  terminations  of  the  a-Al203(0001)  surface  were  examined.  Fig.  6  dis¬ 
plays  a  repeated  view  of  a  unit  cell  with  relaxed  atomic  coordinates  for  a 
Ni(l  1 1)/Al2 03(0001)  interface  calculation.  The  lighter  atoms  in  the  top 
layer  correspond  to  aluminum  and  the  darker  atoms  are  oxygen.  In  the 
Al-terminated  a-Al203(0001)  surfaces,  it  appeal’s  that  interfacial  bonding 
between  the  alumina  and  Ni  substrate  decreases  with  increasing  AI2O3 
thickness.  This  effect  may  contribute  to  the  increased  spallation  observed 
experimentally  with  the  thickening  of  the  TGO. 
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Figure  6  Ni(l  1 1  )/Al203(0001)  relaxed  interface  geometry 


Kruse  et  ft/.194  195  studied  Nb  monolayers  (and  multiple  Nb  layers  to 
simulate  bulk-like  properties)  on  the  O-terminated  <x— Al203(0001)  surface 
and  found  that  Nb  bonded  to  the  hollow  sites  above  the  vacant  octahedral 
sites.  These  authors  emphasized  the  importance  of  ion  relaxations  when 
comparing  competing  adsoiption  geometries,  a  feature  often  lacking  in  many 
reported  studies. 

We  have  also  undertaken  studies  of  the  Ni/ZrCU  interface  again  using 
ultrasoft  pseudopotentials  within  the  VASP  code.178- 196  Up  to  three  layers  of 
cubic  Zr02(l  1 1)  were  “deposited”  onto  the  Ni(l  1 1)  substrate.  As  we  found 
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for  the  Ni/Al203  interface,  the  work  of  separation  (or  adhesion  strength)  is 
quite  large  for  the  first  monolayer;  and,  as  more  Zr02  is  added,  the  intra¬ 
ceramic  bonds  increase  in  strength  at  the  expense  of  the  interfacial  bonding. 
The  adhesion  strength  is  halved  by  the  time  that  3  layers  of  Zr02  arc  present. 
We  find  significant  atomic  distortions  at  the  interface  and  a  suppression 
of  magnetism  in  the  Ni  layers  nearest  to  the  interface.  The  decrease  in 
adhesion  as  Z1O2  grows  on  Ni  is  consistent  with  the  need  for  a  bond  coat 
alloy,  since  our  predictions  suggest  thick  Z1O2  films  may  not  readily  adhere 
to  the  Ni-based  superalloy  under  the  TBC. 

6.2. 2. 2.  Aplications  of  D FT  to  Ceramics  We  have  also  studied  bulk  Z1O2 
in  the  monoclinic,  tetragonal,  and  cubic  phases  as  well  as  their  associated 
low-index  surfaces  within  the  LDA  to  DFT.48  We  found  that  the  most  stable 
surfaces  of  tetragonal(t)  and  monoclinic(m)  zirconia  -  the  t(l  1 1)  and  m(l  11) 
surfaces  -  arc  nearly  degenerate  in  their  surface  energies.  This  provides 
an  explanation  for  the  preferential  stabilization  of  t-ZrC>2  in  small  particles, 
where  surface  energies  rather  than  bulk  cohesive  energies  should  determine 
the  preferred  structure.  By  using  a  Wulff  construction197  and  known  ori¬ 
entation  relationships  between  the  t  and  m  phases,  we  showed  that  t(lll) 
surfaces  will  transform  to  halfm(Ill)  and  halfm(lll)  surfaces  upon  phase 
transformation.  While  the  infill)  surfaces  are  low  in  energy,  the  m(lll) 
surfaces  are  predicted  to  have  somewhat  higher  energies.  Thus,  we  sug¬ 
gested  that  the  suppression  of  the  tetragonal  to  monoclinic  phase  transition 
in  small  particles  may  be  due  to  the  thermodynamically  unfavorable  nature 
of  m(l  1 1)  surfaces  that  would  be  forced  to  form.  This  idea  prompted  us  to 
consider  ways  to  keep  Z1O2  in  the  tetragonal  phase  over  a  very  wide  tem¬ 
perature  range.  One  possible  solution  is  to  embed  the  small  Z1O2  particles 
in  an  alumina  matrix,  for  example,  as  in  the  idea  behind  the  nanolaminate 
films  comprised  of  Zr02  and  AI2O3  multilayers. 

We  then  undertook  a  DFT-GGA  study  (within  the  Projector  Augmented 
Wave  formalism)198  of  the  Zr02-Al203  interface,  in  order  to  understand 
how  alumina  might  serve  this  purpose  of  confining/stabilizing  small  t-Zr02 
particles  and  to  shed  light  on  the  nature  of  the  interaction  between  the  bond 
coat  and  the  TBC  that  occurs  when  the  bond  coat  is  oxidized  to  AI2O3.199 

Our  main  conclusions  are  that  the  alumina/zirconia  interface  is  weakly 
interacting:  we  find  negligible  charge  transfer,  no  evidence  of  covalent 
bonding,  and  a  very  low  adhesion  energy  of  1.065  J/m2.  The  low  adhesion 
energy  is  probably  due  to  the  fact  that  these  surfaces  of  Zr02  and  AI2O3 
reconstruct  to  obtain  approximate  coordinative  saturation,  and  therefore  the 
lack  of  dangling  bonds  on  these  surfaces  minimizes  the  interaction  they  have 
between  them.  This  suggests  that  the  role  of  the  AI2O3  in  the  nanolaminate 
coatings  is  simply  to  act  as  a  physical  barrier  to  growth  of  the  Zr02  layer  and 
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that  there  is  no  true  chemical  bonding  between  these  layers.  Further,  this 
weak  interaction  has  broad  implications  for  thermal  barrier  coatings.  When 
the  bond  coat  oxidizes,  it  is  known  that  AI2O3  forms  and  that  the  lifetime 
of  the  TBC  is  tied  to  the  oxidation  kinetics  of  the  bond  coat.  Therefore 
a  possible  microscopic  explanation  is  now  available:  once  AI2O3  forms, 
there  is  only  a  weak  interaction  between  Zr02-Al203,  and  thus  the  Z1O2 
de-adheres  at  that  point.  A  solution  to  this  problem  may  be  to  minimize 
the  amount  of  A1  in  the  bond  coat.  A1  may  oxidize  more  readily  than  other 
elements  present  in  the  bond  coat,  and  so  perhaps  increasing,  for  example, 
the  amount  of  Cr  in  the  bond  coat  may  help  inhibit  spallation. 

To  summarize  the  application  of  DFT  methods  to  M/C  and  ceramic  inter¬ 
faces,  a  significant  “noise  level”  is  apparent  in  predicted  energetic  features, 
originating  from  different  basis  set  expansions,  different  degrees  of  ionic 
relaxation  allowed  for,  and  different  approximations  for  electron  exchange- 
correlation  effects.  However,  qualitative  trends  from  different  studies  are 
most  often  in  agreement  with  each  other  and  experiments.  We  now  turn  to 
a  more  approximate  DFT  technique  which  will  allow  larger  systems  to  be 
studied. 


6.3.  The  Harris  Functional 


A  conceptually  important  approximation  within  DFT  is  the  Harris  func¬ 
tional.200-201  In  the  Harris  functional,  the  electronic  ground  state  density  p 
in  Eq.  (1 1)  is  approximated  by  a  sum  of  overlapping  atom- like  densities  pm, 
suitably  chosen202-203  for  each  species  m  =  1  •  •  •  M 

M  Nm 

p(r)  =  II  Pm(f  Rmi)  (12) 

m=l  i=1 


where  Rm^  designates  the  position  of  the  itH  atom  of  species  m.  This  den¬ 
sity  ansatz  makes  it  possible  to  express  the  cohesive  energy  and  Hamiltonian 
matrix  elements  in  terms  of  atomic  contributions,  which  is  very  conve¬ 
nient  for  calculational  and  interpretive  purposes.  The  Harris  functional  has 
been  used  as  the  starting  point  for  deriving  many  approximate  total  energy 
schemes.201, 204-206  Inserting  Eq.  (12)  into  Eq.  (11)  we  get 


E[p]=  Y.  £i[v[p]]- I  v[p]pdr  +  Hclassical  [p]  +  Uxc[p]  +  Tion  (13) 

iEoccup 

where  the  one-electron  Schrodinger  equations  in  the  Kohn-Sham  formal¬ 
ism161  was  used  to  get  the  first  two  terms  and  the  one-electron  potential  v  is 
given  by 

6(UXC  +  Uciass-LCQ[J 

6p  lp  =  p 


v  p  = 


(14) 
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The  main  problem  with  the  Harris  functional  is  that  charge  transfer  and 
screening  effects  arc  not  included  properly,  due  to  the  static  density  ansatz, 
Eq.  (12).  This  has  led  Finnis,  on  the  basis  of  the  image  charge  model207,208 
for  M/C  bonding,  to  object  that  the  Harris  functional  -  and  therefore  ap¬ 
proximate  total  energy  schemes  derived  from  the  Harris  functional,  such  as 
non-self-consistent  tight  binding  -  are  unsuitable  for  describing  M/C  bond¬ 
ing.  He  asserts  that  the  attractive  component  of  the  M/C  cohesion  is  mainly 
a  polarization  effect  in  the  image  charge  model,  which  we  will  discuss  later. 
However,  despite  the  static  density  ansatz,  Eq.  (12),  static  screening  effects 
are  included  via  the  eigenvalue  sum  in  Eq.  (13).  Thus  it  remains  to  be 
demonstrated  on  a  broad  range  of  M/C  systems  how  well  the  Harris  func¬ 
tional  performs.  Encouraging  results  have  been  reported  for  Ag/MgO  and 
Al/MgO  with  and  without  C  and  S  impurities.177,209 

6.4.  Tight  Binding 

Tight  binding  (TB)  schemes210,211  refer  to  a  rather  ill-demarcated  group 
of  electronic  structure  techniques  for  describing  interatomic  interactions. 
The  key  merit  of  TB  methods,  compared  to  lower  levels  of  interatomic 
potentials,  is  that  they  include  quantum  effects  explicitly,  although  at  a 
fairly  approximate  level.  Their  simpler  nature  produces  a  large  computa¬ 
tional  speed  gain  over  ah  initio  methods.  This  enables  approximate  studies 
of  electronic  structure  effects  in  extended  systems  such  as  TBC’s.  In  the 
chemically-oriented  literature,  TB  methods  are  associated  with  Extended 
Hiickel  Theory(EHT)212,213  or  variants  like  Atom  Superposition  and  Elec¬ 
tron  Delocalization  Molecular  Orbital  (ASED-MO)  theory,  which  is  EHT 
augmented  by  an  inter-nuclear  repulsion  term.214 

The  archetype  for  the  TB  total  energy  expression,  ETB,  for  ionic  coordi¬ 
nates  {Rn|  is 

ETB{Rn}=  Y-  £i{Ru}  +  Erep{Ru}  (15) 

iGoccup 

where  Erep  is  a  repulsive  term  balancing  the  propensity  to  gain  energy  in  the 
spectral  sum  by  increasing  orbital  coupling.  The  electronic  spectrum  {e-J  is 
obtained  by  solving  a  matrix  eigenvalue  problem 

|H  -  SE\  =  0  (16) 

corresponding  to  an  implicit  or  explicit  atomic  basis,  usually  chosen  to  be 
minimal,  i.e.  one  s-orbital,  three  p-orbitals  per  main  group  atom  and  five 
d-orbitals  per  atom  for  transition  metals.  H  and  S  arc  the  Hamiltonian  and 
overlap  matrices  corresponding  to  this  basis.  Comparing  Eqs.  (15)  and  (16) 
to  Eq.  (13),  the  inherent  approximations  of  TB  appeal-  :  (i)  The  overlap 
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matrix  S  is  often  set  to  the  unit  matrix,  i.e.  the  atomic  basis  is  assumed 
orthogonal  (orthogonal  TB)  -  this  reduces  accuracy.  Proper  inclusion  of  the 
overlap  matrix  S  (non-orthogonal  TB)  increases  transferability,  but  is  com¬ 
putationally  more  demanding,  (ii)  The  Hamiltonian  and  overlap  matrices 
arc  often  fitted  to  simple  pairwise  analytic  forms  from  a  limited  data-set, 
instead  of  being  appropriate  matrix  elements  with  respect  to  a  given  basis. 
This  reduces  accuracy,  (iii)  The  repulsive  term  ETep  is  usually  represented 
by  an  isotropic  pair-potential,  summed  over  all  atomic  pairs.  The  coun¬ 
terpart  for  ETep  in  Eq.  (13)  is  nonlinear  in  the  exchange-correlation  terms, 
therefore  accurate  representation  of  ETep  requires  environment  dependent 
pair-potentials215-217  or  three-body  corrections,  (iv)  Usually  there  is  a  lack 
of  charge  self-consistency  (i.e.,  charge  redistribution,  which  is  neglected 
in  the  density  ansatz,  Eq.  (12)).  This  is  usually  only  done  by  adjusting 
the  diagonal  Hamiltonian  matrix  elements  Hy,.211,218  Charge  transfer  may 
induce  long-ranged  effects,  i.e.  they  may  create  a  Madelung  potential  or 
induce  band-bending.  The  most  modern  TB  schemes  put  emphasis  on  a 
more  realistic  and  detailed  inclusion  of  self-consistency  effects.206,219 

The  accuracy  to  which  these  points  are  treated  determines  the  transferabil¬ 
ity  of  the  TB  scheme.  For  low  level  TB  schemes,  at  most  only  the  qualitative 
trends  should  be  trusted.  A  pitfall  in  this  context  is  overly  parameterized 
TB  schemes,  i.e.  too  many  fitted  input  parameters  -  this  strategy  will  tend  to 
hide  inherent  limitations.  In  principle,  care  fully  worked  out  TB  schemes  can 
achieve  almost  the  accuracy  of  self-consistent  LCAO-DFT204, 206  without  use 
of  empirical  fitting  parameters. 

Many  interesting  TB  studies  have  been  performed  on  the  energetics  and 
electronic  structure  of  M/C  interfaces.  Most  TB  studies  have  focused  on 
transition  metals  deposited  on  a— alumina.  Alemany  et  al. 220  obtained  band 
structures  using  TB  in  an  extended  Hiickel  framework  to  study  a  variety  of 
crystalline  alumina  terminations,  where  O  anions  were  in  the  O2-  charge 
state.  They  generally  found  attractive  interactions  between  adsorbed  3d- 
transition  metal  atoms  coordinating  with  Al3  1  cations  but  repulsive  interac¬ 
tions  when  coordinating  with  O2-  anions,  where  the  repulsiveness  increased 
with  nobleness.  They  found  only  small  steric  effects  in  this  pattern,  which 
provides  support  for  phenomenological  approaches  to  quantifying  interface 
cohesion.  They  attributed  the  transition  metal-Al3+  attraction  to  mixing 
of  the  adsorbate  d-states  and  the  coordinatively  unsaturated  Al3+  dangling 
bond  (empty  sp-hybrid)  surface  states. 

When  O  ions  were  in  the  O-  charge  state,  the  anions  were  found  to 
be  very  reactive,  with  maximum  interface  adhesion  for  V  and  decreasing 
monotonically  towards  each  side  in  the  d  series,  but  with  Cu  anomalously 
stabilized.  The  interface  strength  of  0“  terminated  a— Al203(0001)  is  esti- 
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mated  to  twice  that  of  the  Al3+  terminated.  This  picture  is  in  some  qualitative 
disagreement  with  the  findings  of  Nath  and  Anderson,221  who  found  the  in¬ 
terface  bond  strength  to  0“-terminated  <x— Al203(0001)  decreasing  mono- 
tonically  with  nobleness  along  the  entire  transition  series.  The  Nath  and 
Anderson  study  used  ASED-MO  theory  to  investigate  their  semiempirical 
cluster  models  where  bonding  was  composed  of  covalent,  charge-transfer, 
and  two-body  repulsion  terms.  The  uncertainty  between  different  studies 
even  at  the  qualitative  level  is  a  signature  of  the  sensitivity  on  TB  input 
parameters  and  specific  approximations  in  the  TB  scheme  applied.  This  is 
also  seen  comparing  the  results  of  Anderson  et  al.  ’s222  ASED-MO  cluster 
calculations  to  Ward  et  al.  ’s223  approximate  weighted  Hq  extended  Hiickel 
slab  calculations  for  Pt  on  a-Al203(0001).  Both  authors  agree  that  Pt  binds 
to  the  O2- -terminated  surface  by  ~  0.1  eV/Pt,  but  for  the  Al3+ -terminated 
surface  Ward  et  al.223  find  3.7  eV/Pt  whereas  Anderson  et  al  222  obtain  2.5 
eV/Pt.  Interestingly,  in  the  same  study.  Ward  et  al.  find  that  Rh  does  not 
bind  at  all  whereas  Pd  does,  opposite  to  the  nobleness  trend  found  for  the 
O-  termination. 

Ohuchi  and  Kohyama191  conducted  empirical  TB  energy-band  (slab) 
calculations  for  many  transition  metal/alumina  interfaces,  where  the  a- 
Al203(0001)  was  Al3+  terminated.  They  generally  found  the  interface 
cohesion  to  change  from  ionic  to  covalent  character,  when  moving  from 
left  to  right  in  the  transition  series.  They  focused  especially  on  Nb/a- 
Al203(0001)  interfaces,  where  they  noticed  that  the  electronic  perturbation 
was  essentially  confined  to  the  Nb-layer  nearest  to  the  AI2O3  surface.  This 
observation  is  in  contrast  to  the  DFT  study  by  Kruse  et  al.,195  which  found 
that  the  electronic  perturbation  of  the  Nb  layers  was  significant  away  from 
the  nearest  layer.  Although  the  AI2O3  in  the  study  by  Kruse  et  al.  was 
O-terminated,  this  is  more  likely  to  be  a  consequence  of  truncation  of  the 
TB  Hamiltonian  to  nearest-neighbor  contributions  and  limited  inclusion  of 
self-consistency  in  the  work  of  Ohuchi  and  Kohyama. 

Impurities  and  dopants  at  the  M/C  interface  arc  very  important  in  TBC 
systems.  This  issue  also  has  been  addressed  by  TB  studies.  For  the  Ni/Al203 
system,  Hong  et  al.224  using  clusters  treated  with  ASED-MO  theory,  found 
S  impurities  significantly  decrease  the  interface  strength,  due  to  S_2-0-2 
closed-shell  repulsions  increasing  the  distance  of  the  Al3+  from  the  Ni 
surface  layer.  Anderson  et  al.225  found  increased  interface  adhesion  for  Y 
dopants  at  the  Ni/Al203  interface.  Open  d-shell  transition  metals  are  indeed 
expected  to  have  the  greatest  bonding  to  ceramics.  Further,  Ni  and  Y  form 
rather  stable  alloys,  suggesting  favorable  Ni-Y  interaction  at  the  interface. 
Application  of  dopants  for  TBC’s  of  course  requires  the  dopants  to  be  stable 
at  the  interface  upon  thermal  cycling. 
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The  moral  of  the  TB  story  is  to  judge  the  results  with  some  caution,  since 
as  we  outlined  above,  even  predicted  qualitative  trends  vary,  depending  on 
the  parametrization  and  particular  form  of  the  TB  scheme  applied. 

6.5.  Semiempirical  Potentials 

Many  issues  important  for  the  understanding  and  engineering  of  real  life 
interfaces  involve  long  length  scales,  such  as  static  and  dynamic  macro¬ 
scopic  properties,  interface  defects,  toughness,  wear,  corrosion  resistance, 
complex  reaction  phases,  phase  transitions,  and  interface  deadhesion  pro¬ 
cesses.  These  aspects  of  interface  modeling  are  far  outside  the  reach  of 
ab  initio  and  TB  methods  described  above,  which  currently  -  and  for  years 
to  come  -  are  limited  to  studying  model  systems  mimicking  authentic  M/C 
interfaces. 

Therefore,  an  important  objective  in  the  theory  branch  of  M/C  interface 
research  is  to  envisage  new  types  of  interatomic  potentials  suitable  for  mod¬ 
eling  M/C  interfaces  to  bridge  the  length  scale  gap  between  theory  and 
experiments.  In  this  context,  ab  initio  studies  are  very  useful:  first  for  ex¬ 
ploring  the  qualitative  nature  of  M/C  chemistry  which  must  by  captured  by 
novel  model  interaction  potentials;  secondly,  to  provide  a  broad  database  for 
testing  new  interaction  potentials. 

For  metals  and  ceramics  considered  separately,  reliable  interaction  poten¬ 
tials  have  been  developed  in  the  past.  For  many  metals,  embedded  atom  type 
potentials226  230  have  proven  successful.  Alternatives  and  refinements  have 
been  developed.231-233  For  ceramics,  rigid-ion  and  shell  model  potentials234 
and  their  refinements235-238  have  proven  equally  capable. 

For  the  M/C  interface,  atomistic  interaction  potentials  have  developed 
around  the  idea  of  the  image  charge  interaction,239  where  the  ceramic  ions 
arc  attracted  to  their  electrostatic  images  induced  in  the  metal  (assumed 
perfectly  conducting).  This  simple  picture  gives  the  right  order  of  magnitude 
for  the  cohesion.  Image  charge  interaction  is  only  accurate  beyond  the 
distance  where  wavefunctions  of  metal  and  ceramic  overlap.  In  that  case,  the 
ceramic  acts  as  an  external  Coulombic  potential  on  the  metal  and  Finnis208 
derived  the  image  charge  interaction  from  DFT  for  this  case.  At  close 
range,  where  wavefunctions  overlap,  orbital  mixing  contributes  to  the  M/C 
interaction.  Finnis208  argues  that  this  may  be  represented  by  a  closed-shell 
repulsion  term  by  considering  the  ceramic  as  a  perturbed  noble  gas  crystal, 
where  protons  have  been  transferred  from  anions  to  cations  but  the  electron 
distribution  is  kept  frozen.  This  point  of  view  may  be  valid  for  weak  M/C 
interfaces,  but  needs  corrections  for  strongly  interacting  M/C  interfaces  with 
significant  charge  transfer  or  covalent  contributions  to  the  bonding. 
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An  interesting  implication  of  the  image  interaction  picture  is  that  metals 
will  generally  adhere  more  strongly  to  polar  rather  than  nonpolar  surfaces. 
This  is  because  the  net-charged  ceramic  surface  will  induce  a  net-charged 
image  in  the  metal.  By  contrast,  polar  ceramic  surfaces  themselves  are 
unstable,  also  for  electrostatic  reasons. 

A  refinement  of  the  image  model,  the  discrete  classical  model  (DCM), 
was  also  formulated  by  Finnis.208  Here  the  perfectly  conducting  halfspace 
representing  the  metal  is  replaced  by  a  halfspace  of  conducting  spheres.  This 
discretization  of  the  metal  is  necessary  for  atomistic  modeling;  and  further 
this  refinement  gives  a  more  realistic  short-range  description  than  the  image 
model,  where  the  interaction  energy  diverges  at  the  image  plane  (i.e.,  the 
metal  surface).  Unfortunately,  Duffy  et  al.240  found  that  the  DCM  predicts 
that  Ag  adsorbs  over  Mg+2  instead  of  O-2  on  the  MgO(OOl)  surface,240  at 
variance  with  all  electronic  structure  calculations. 

A  hybrid  model  for  full  atomistic  relaxation  was  used  recently  by  Purton 
et  al.241  to  simulate  the  weakly  interacting  Ag/MgO(001)  interface.  The 
metal  was  described  by  an  embedded-atom  potential,  the  ceramic  by  a  shell 
model  and  the  M/C  interaction  by  the  DCM  (i.e.  an  atom  interacts  via  2 
different  potentials).  These  authors  found  that  Ag  adsorbed  over  O-2  and 
attributed  this  to  the  fact  that  they  used  a  different  short-range  potential  in 
the  DCM  than  Duffy  et  al.240  used. 

This  illustrates  that  we  arc  far  from  the  point  where  we  can  develop  robust 
and  transferable  semiempirical  potentials  for  the  general  M/C  interface.  This 
task  is  fundamentally  difficult  because  the  bonding  changes  character  at  the 
interface,  from  metallic  to  ionic,  and  because  of  the  structural  complexity 
present.  This  probably  requires  explicit  inclusion  of  quantum  effects,  at  least 
at  a  low  level. 

Because  semiempirical  potentials  for  M/C  interfaces  arc  not  yet  a  mature 
field,  simulation  results  that  rely  critically  on  details  of  the  metal  and  ceramic 
layer  adjacent  to  the  interface  need  to  be  judged  with  great  caution,  because 
of  the  limited  transferability  expected  for  these  potentials.  For  strongly 
coupled  M/C  interfaces,  no  reliable  functional  forms  of  atomistic  interaction 
potentials  have  been  proposed  yet. 

7.  MULTISCALE  MODELING 

As  mentioned  before,  a  salient  problem  with  the  atomistic  modeling  of 
realistic  M/C  interfaces  is  the  contemporary  presence  of  effects  on  multiple 
length  scales,  chemical  interactions  in  the  interface  region  along  with  de¬ 
fect  structures  and  strain  fields  extending  far  away  from  the  interface  region. 
Naturally,  this  is  a  well-recognized  difficulty  in  many  modeling  applications, 
and  numerous  techniques  have  been  introduced  that  attempt  to  overcome  this 
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problem.  One  useful  approach  is  to  couple  an  active  region  (in  our  case  the 
interface  vicinity)  which  is  treated  explicitly  atomistically,  together  with  the 
surrounding  system  within  linear  elasticity  theory.242,243  Course-Grained 
Molecular  Dynamics  (CGMD)  is  also  designed  to  allow  the  treatment  of 
larger  systems  than  would  be  accessible  with  conventional  ab  initio  tech¬ 
niques.244  In  CGMD,  a  statistical  course-graining  method  is  used  to  derive 
the  equations  of  motion  from  finite  temperature  Molecular-Dynamics(MD). 
Accordingly,  these  equations  coincide  with  conventional  MD  for  atomic- 
scale  mesh  sizes,  which  allows  a  smooth  coupling  of  the  length  scales.  As 
a  result,  CGMD  provides  a  means  to  computationally  concentrate  on  “ac¬ 
tive  regions,”  by  creating  finer  meshes  there,  while  permitting  large-system 
simulations. 

Another  promising  approach  is  the  quasicontinuum  method  proposed  by 
Tadmor  el  a l.,245  which  combines  the  finite  element  method  with  adaptive 
mesh  techniques.  Atoms  arc  partitioned  into  “local”  and  “nonlocal”  atoms, 
according  to  the  strain  at  the  atom.  Atoms  characterized  by  low  strain  (local 
atoms)  are  typically  far  from  the  active  region,  and  the  energy  of  these  atoms 
is  well-approximatedby  atoms  in  a  perfect  crystal  subjected  to  a  similar  (low) 
strain.  This  energy  can  be  obtained  from  standard  first  principles  methods. 
Technically,  the  energy  of  the  local  atoms  in  the  large,  remote  regions  is 
obtained  by  finite  element  integration  on  a  sparse  grid.  The  nonlocal  atoms 
are  subjected  to  high  strain  in  the  vicinity  of  an  active  region.  The  cohesive 
energy  of  a  “nonlocal”  atom  is  calculated  as  an  explicit  function  of  its 
coordination  shell  in  real  space,  out  to  a  given  cutoff  radius,  which  is  more 
time-consuming.  Conventional  atomistic  simulations  normally  consider  all 
atoms  as  “nonlocal”,  which  is  why  the  quasicontinuum  method  can  treat 
of  order  1000  times  more  atoms  than  feasible  with  conventional  atomistic 
simulation. 

A  limitation  of  the  present  formulation  is  that  it  relies  on  the  partitioning 
of  the  energy  into  atomic  contributions,  i.e. 

All  atoms 

Etot  =  Y_  El  (IV) 

i=i 

where  E-r  is  an  average  energy  for  a  “representative”  atom.  This  is  a  virtue 
of  approximate  schemes  only,  like  the  embedded  atom  method  and  pair- 
potentials;  the  exact  total  energy  is  a  collective  quantity  and  imposing  energy 
separability,  Eq.  (17),  for  first  principles  methods  introduces  some  degree  of 
ambiguity,  especially  for  the  atoms  treated  as  “nonlocal”.  Also  due  to  the 
energy  separability  requirement,  extension  of  the  quasicontinuum  method  to 
multicomponent  systems  is  less  clear,  if  energies  and  forces  are  generated 
from  ab  initio  methods.  One  may  be  able  to  get  around  this  problem  by 
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defining  E-l  in  the  sum  in  Eq.(17)  to  be  the  total  energy  per  unit  cell  for  an 
associated  strained  reference  system. 

No  results  for  M/C  interfaces  using  these  rather  new  multiscale  modeling 
techniques,  even  with  approximate  total  energy  schemes,  have  emerged  yet 
to  our  knowledge;  but  studies  of  weakly  coupled  M/C  interfaces  using  a 
hybrid  potential,  like  that  of  Purton  et  al .,241  would  certainly  be  feasible. 

In  addition  to  bridging  length  scales,  as  in  the  quasicontinuum  method,  it 
is  equally  important  to  consider  how  to  bridge  time  scales.  For  example,  the 
time  scale  of  atomic  motion  is  many  orders  of  magnitude  smaller  than  the 
time  scale  for  materials  creep.  One  route  to  accomplish  time  scale  bridging 
is  the  Kinetic  Monte  Carlo  (KMC)  technique.246  In  this  method,  the  system 
is  evolved  stochastically  according  to  probabilities  given  by  ratios  of  rate 
constants  lq  to  kref ,  where  kTef  is  the  fastest  rate  constant  in  the  system.  By 
first  choosing  randomly  a  species  in  the  system  and  then  choosing  randomly 
from  a  palette  of  possible  elementary  events  that  may  change  the  nature  of 
that  species,  and  comparing  the  rate  constant  ratio  for  the  selected  elementary 
event  to  a  random  number  between  0  and  1,  one  has  a  prescription  for  how 
to  evolve  the  system  in  a  kinetic  rather  than  an  equilibrium  manner.  While 
this  provides  the  means  to  get  to  long  time  scales,  it  does  have  one  major 
drawback.  One  must  specify  in  advance  the  allowed  elementary  processes; 
in  so  doing,  one  may  actually  constrain  the  behavior  of  the  material  under 
study.  Further,  deterministic  dynamics  have  been  sacrificed  at  the  altar  of 
stochastic  behavior  in  order  to  gain  the  increase  in  time  scales.  However,  first 
principles-derived  KMC  has  been  used  successfully  to  study  kinetic  behavior 
for  a  number  of  different  materials  by  our  group247  and  others248-252  (not  all 
of  which  arc  listed  here). 

Recently,  Voter  proposed  a  method  known  as  hyperdynamics,  or  hyper- 
MD,  designed  to  accelerate  molecular-dynamics  simulations.253'254  This 
acceleration  is  achieved  by  including  a  bias  potential  that  raises  the  potential 
energy  in  the  valley  regions  of  the  potential  energy  surface  but  leaves  the 
potential  in  transition  state  regions  unaffected.  Because  the  bias  potential 
does  not  affect  transition  state  regions,  even  though  the  Transition  State 
Theory  (TST)  escape  rate  is  enhanced,  the  ratio  of  TST  escape  rates  from 
a  given  state  to  adjacent  states  is  preserved.  As  a  result,  the  system  should 
advance  sequentially  at  an  accelerated  pace  while  preserving  the  relative 
probabilities  of  exact  dynamics.  Provided  that  the  assumptions  inherent  in 
TST  are  obeyed  for  the  processes  of  interest  (i.e.,  no  saddle  point  recrossings 
and  no  correlated  events),  this  is  an  attractive  method  for  reaching  time  scales 
orders  of  magnitude  longer  than  those  achieved  by  conventional  molecular 
dynamics.  Voter  has  also  proposed  a  method  to  effectively  parallelize  MD 
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simulations.255  This  technique  and  the  hyper-MD  arc  beneficial  in  decreasing 
simulation  time  necessary  for  the  observation  of  infrequent  events. 

Kinetic  Monte  Carlo  and  hyperdynamics  methods  have  yet  to  be  applied  to 
processes  involved  in  thermal  barrier  coating  failure  or  even  simpler  model 
metal-ceramic  or  ceramic-ceramic  interface  degradation  as  a  function  of 
time.  A  hindrance  to  their  application  is  lack  of  a  clear  consensus  on  how 
to  describe  the  interatomic  interactions  by  an  analytic  potential  function. 
If  instead,  for  lack  of  an  analytic  potential,  one  must  resort  to  full-blown 
density  functional  theory  to  calculate  the  interatomic  forces,  this  will  become 
the  bottleneck  that  will  limit  the  size  and  complexity  of  systems  one  may 
examine,  even  with  multiscale  methods. 

8.  PHENOMENOLOGICAL  APPROACHES 

Phenomenological  approaches  has  been  very  successful  in  some  areas, 
e.g.  Miedema  theory256  for  predicting  many  quantities  in  metallurgy.  The 
essential  task  in  phenomenological  theories  is  to  identify  a  suitable  set  of 
physically  meaningful  variables,  which  are  “linearly  independent”,  to  char¬ 
acterize  the  materials.  Experimental  data  are  then  correlated  against  this  set 
of  variables  and  functional  relations  are  fitted. 

Li257  explored  such  possibilities  for  M/C  interfaces.  He  found  that  the 
interface  cohesion  correlated  well  with  average  metal  electron  density  in 
bond  regions,  the  free  energy  of  formation,  the  band  gap  and  the  conductivity 
of  the  ceramic.  This  study  mainly  focused  on  transition  metals,  but  the 
correlation  with  electron  density  seemed  also  to  encompass  simple  metals 
and  semiconductors  as  interface  partners.  Chatain  el  n/.258  stressed  the  local 
chemical  aspect  by  correlating  interface  cohesion  with  a  linear  combination 
of  the  solution  enthalpy  of  the  ceramic  anion  and  cation,  respectively,  in  the 
bulk  metal,  plus  a  van  der  Waals  term.  For  simple  M/C  interfaces,  Bordier 
and  Noguera’s  model  tight  binding  study259  indicates  that  the  Fermi  level 
of  the  metal  and  the  degree  of  ionicity  of  the  ceramic  should  be  important 
variables  for  correlating  experimental  data.  Of  course,  the  diversity  of  the 
approaches  and  variable  sets  above  illustrates  that  there  is  no  unique  strategy 
for  understanding  M/C  cohesion. 

9.  CONCLUSIONS  AND  OUTLOOK 

Understanding  the  behavior  of  the  interfaces  and  bulk  materials  involved 
in  thermal  barrier  coating  failure  due  to  the  extreme  environment  created  in 
aircraft  engines  is  still  in  its  infancy.  This  is  primarily  because  the  system 
involves  complex  interfacial  chemistry  and  the  materials  issues  span  large 
length  and  time  scales.  In  this  review,  we  have  focused  on  the  atomic  level 


40  Asbjom  Christensen,  Emily  A.  Asche  and  Emily  A.  Carter 


characterization.  Once  that  is  specified,  it  will  be  imperative  to  draw  links 
between  the  atomic  and  the  microstructural  scales  in  order  to  understand  the 
materials  failure  mechanisms  completely. 

Ultimately,  an  important  goal  of  the  structural  studies  described  above  is 
to  provide  a  foundation  for  experiment  and  theory  to  investigate  the  dynamic 
behavior  of  M/C  interfaces,  especially  under  conditions  similar  in  severity 
to  those  found  in  the  combustion  region  of  a  gas  turbine  engine.  However,  as 
mentioned  above,  theoretical  dynamical  studies  of  this  type  would  be  highly 
nontrivial  and  to  date  arc  nonexistent.  In  fact,  constructing  techniques  to 
make  these  dynamical  studies  feasible  constitutes  an  active  area  of  our  own 
research  program.  Why  have  such  simulations  not  been  performed?  It  is 
because  of  the  -  at  best  -  murky  understanding  of  M/C  bonding.  In  other 
words,  the  lack  of  dynamics  simulations  on  such  systems  is  indicative  of  a 
much  deeper,  unanswered  question;  namely,  what  is  the  true  nature  of  the 
interaction  between  metal  atoms  and  the  atoms  present  in  a  ceramic?  It  is  a 
much  more  difficult  and  complex  system  to  characterize  than  either  a  metal 
itself  or  an  isolated  ceramic,  each  of  which  can  be  described  independently 
with  much  simpler  models  (e.g.,  glue  or  embedded-atom  models  for  metals, 
ionic  or  shell  models  for  ceramics,  etc.).  It  is  the  interface  interaction  which 
eludes  understanding  to  date.  Thus,  the  emphasis  in  this  review  on  static 
rather  than  dynamic  properties  of  these  materials  points  out  that  theorists 
arc  not  yet  in  a  position  to  examine  these  systems  dynamically  on  the  large 
length  and  time  scales  which  may  prove  necessary  to  fully  simulate  thermal 
barrier  coating  failure. 

That  said,  we  can  still  glean  some  insight  from  the  static  investigations, 
enough  to  even  make  some  tentative  conclusions  about  how  to  improve 
TBC’s.  It  is  clear  that  the  ideal  coating  should  consist  of  a  stress-tolerant 
ceramic  and  an  oxidation-resistant  bond  coating  with  nearly  the  same  co¬ 
efficient  of  thermal  expansion  and  elastic  modulus  as  both  the  ceramic  and 
the  Ni-based  superalloy.  The  cost  savings  that  such  an  ideal  coating  would 
produce  is  nontrivial:  it  has  been  estimated  at  107  gallons  of  fuel  per  year 
for  a  fleet  of  250  airplanes. 

Use  of  nanostructured  pure  Z1O2  might  help  create  such  an  ideal  coating, 
as  recent  experiments  have  shown  that  one  can  stabilize  the  tetragonal  phase 
of  undoped  Z1O2  up  to  a  critical  grain  radius  of  6  nm.49  We  have  provided 
an  explanation  for  the  special  stability  of  t-ZrCU  in  nanocrystals,  based  on 
surface  energy  arguments.  Thus,  appropriately  synthesized  nanocrystals 
of  Zr02  might  allow  the  use  of  pure  rather  than  doped  zirconia,  which 
would  eliminate  possible  sources  of  materials  failure  due  to  dopant  diffusion. 
Furthermore,  use  of  t-ZrCU  would  provide  the  stress  tolerance  required  via 
the  transformation  toughening  mechanism  observed  for  t-ZrCU-48 
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It  may  be  that  the  life  expectancy  of  a  TBC  is  even  more  sensitive  to 
the  bond  coating  composition  than  to  the  ceramic  top  coat  composition.260 
Oxidation  of  the  bond  coating  to  AI2O3  may  also  be  a  significant  source 
of  TBC  spallation.  We  suggest  that  reducing  the  amount  of  A1  in  the  bond 
coat  may  be  necessary  to  avoid  spallation,  based  on  our  first  principles 
calculations. 

Based  on  the  theoretical  work  done  to  date,  it  appeal's  that  rare  earth  or 
transition  metals  with  open  d-  or  f-shells  (the  less  noble  the  better)  show  the 
greatest  adhesion  to  ceramics,  presumably  because  there  can  be  a  degree  of 
covalent  bonding  at  the  interface.  On  the  other  hand,  one  wants  to  utilize 
transition  metals  that  are  either  oxidation-resistant  or  form  oxides  that  adhere 
strongly  to  Zr02-  This  may  explain  the  presence  of  Y  in  the  bond  coating, 
since  oxidation  of  the  Y  metal  to  Y2O3  will  produce  an  oxide  that  forms 
solid  solutions  with  Z1O2,  i.e.  it  is  strongly  interacting.  It  suggests  that  one 
should  consider  metals  that  form  other  cubic-based  metal  oxides,  such  as 
cerium,  for  addition  to  the  bond  coat. 

As  one  can  see,  exploring  the  atomic-level  properties  of  thermal  barrier 
coatings  involves  a  complex  set  of  issues.  We  have  attempted  to  provide 
some  background  on  methods  available  for  preparing  and  characterizing 
TBC’s  experimentally,  as  well  as  presenting  the  hierarchy  of  theoretical 
methods  available  in  computational  materials  science,  including  the  trade¬ 
offs  involved  in  the  use  of  each  one.  It  is  our  hope  that  a  microscopic 
understanding,  provided  both  by  ever-improving  theory  and  experiment,  of 
the  basic  mechanisms  for  metal-ceramic  adhesion,  ion  and  metal  atom  diffu¬ 
sion,  and  ceramic  phase  transformations  eventually  will  lead  to  substantive 
improvements  in  macroscopic  properties  of  TBC’s,  resulting  in  extended 
survival  in  the  extreme  environment  present  in  an  aircraft  engine. 
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